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"STATIC AND DYNAMIC ANALYSIS OF ‘GUY CABLES 


-more exact manner than has been previously accomplished without causing ad-- 


presented in paper was completed in 1955 ‘in an ‘unpublished 
_ study? of flexible towers. It has been available to many engineers | and oo 
- Comments by several engineers resulted in the preparation of this | paper. 
_ though the presentation is directed toward guyed towers, it is believed ae 4 
- formulas will prove useful for other needs such as the analysis of of transmission 7 
lines and unloaded back stays for suspension bridges. 
_ In the design and construction of many structures, suchas towers to support _ 
antennas that radiate information at line-of-sight frequencies, the lateral sup- . 
port is 5 furnished by guys composed of Steel strand cables. _ Because these ca- 


is useful to have accurate and reliable formulas which relate the variables as- _ 
sociated with hanging c: cables, ‘such as the change in cable tension as a function — 


Note.—Discussion open until June extend the closing date one 
7 written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Structural Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No, ST 1, January, 1961 
I prof, and Chmn, of Civ. Engrg., Univ. of Delaware, Newark, Del. —— a 
244 Static and Dynamic Analysis of Tall Flexible Towers,” ” by D. L. Dean, Ph, 
Dissertation, Univ, of 1955. 
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2% a _ Guy cables are analyzed for static and dynamic conditions, and examples — ae 
ive illustrate he formulas derived a sag is aceon d_ for in a 
ae 
j 
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“of the tower movement. In some special cases, formulas developed by consid- 
ering cables as” 3 straight two- -force members are suf sufficient. However, it is u- 
-sually necessary to consider their r sag. 
Tn existing analyses3,4,5,6 of guy cables, the geometry is assumed to be | 
-_ “parabolic. This assumption introduces errors that are quite small in some ] 
instances, such as the formula for cable length if the chord is nearly horizon- 
tal If the chord is not horizontal, however, symmetry (with respect to the ¢ 
chord mid-point) is lost and the approximation introduces significant errors. 
The parabola- -for r-catenary approximation is not made in this paper. _ The for-_ 
_ mulas derived by adhering to the properties of the catenary are as simple as 
those based upon the properties | - the parabola. Therefore, the inaccuracies" 
in the latter are not warranted. + 
In the static analysis that follows, cable movements and load changes are. 
: to occur slowly (as compared with natural frequencies) . In the dy- 
namic analysis, the inertia a of the cables is considered. 
| When the tower vibrates normal to its a axis, the guys experience longitudinal - 
_ and transverse components of the movement. If the tower movement is not in 
- the plane of the guy, there are two components of transverse cable vibration. 
Preliminary computations indicate that the natural frequencies of towers and 
their guys are of the same order of magnitudes for lateral vibration. - There- ; 
fore, in computing the lateral resistance that guys offer to tower movements, 
the inertia of the guys should be considered. The frequencies of the modes of 
- longitudinal vibration of guys are not of the same order of magnitude as those Z 
of towers, and need not be included in the computation of the longitudinal | re-— a 


ratio of the fundamental frequency of longitudinal vibration to lateral 


ration of a is 2 is —, where where i the cable strain. For cables” used as 


€ 
"guys, the is 40. “Thus, t the is 
pared with the fundamental frequencies of ae, so that a static analysis — 7 
is sufficient for longitudinal motion. RY 
te Notation.—The letter symbols adopted for use inthis paper ‘are edefinedwhere _ 
they first appear, in the illustrations or in the text, and are 


ie we Equations and Properties. —The cable is considered to be outer 
4 flexible so that the differential equation of the cable curve can be written by 
summing | the vertical forces which act on an element. Some of the details of 


the derivation of the equation for the cable curve follow. = et 


3 «structural Frameworks,” by | c. tT. Morris, John Wiley an and Sons, ‘New York, Chap- 
= 4 «Length, Tension, and Sag of Stay Ropes; ‘ by zs Walmsley, Inst of dil ey ., No. 
"5 “Design of Multi-Level Guyed Towers,” by E. Cohen and H, Perrin, . Proceedings, 7 
ASCE, Vol. 83, No. ST 5, September, 1957, 
«Amplification of Stress and in R, Rowe, Pro-— 
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GUY CABLES 
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integration constants aj and ag are 


= Oatx=Oandy=hatx=1l 


i 


(4) 


| 

dy 1+ (2) | al 
— luated from tl 

— 


e te msion, T, any x, is F From Eq 


Pt of the other cable variables H, 1, h, and q. “Such an expression will ire 
nish the key relations needed to express changes in one variable 7 terms of 


* 


a 
| [h cot coth 


‘The of cable at tan ‘needed one the computation of 


= sinh [2r+ay] = 
‘The sag, denoteds and measured to one chord ( Fig. 
1), may be desired. Occasionally this quantity is used as a criterion for es- 
_ tablishing the initial erection stress for the guys, ‘such as, s <0.015 L for a 


e. 1/2 in. coating of ice. At Xg, the horizontal distance to the point of a 


dy _ an 
sag, 
dx 


a cosh sin ho - cosh - cosh ay) 
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cabl 
— 
fet e The slope of the cable curve at its upper and lower extremities is frequently eee 
be usedwitha sight- 


Xs sin y - xs) 


dimensionless: r= = , and functions im involving: 


“the ratio ‘between the vertical force necessary to support one- ~half the weight — -@ 


the cable. In guy cables, or cables functioning as guys, r is usually small, 
rarely exceeding 0.5. Therefore, sinh r, cosh r, and tanh r can be replaced = 
a few terms of their respective series without a meaningful loss of accuracy. 
It will often be convenient and helpful to do this in the developments to follow, 7 
and many expressions will be simplified considerably by making such a sub- | 
stitution. An example is the constant a1, V whichcan be rewritten in a form that — 


2) 


Z Eq. 9 can also be changed to forms that make computations easier. This a 
equation, relating the cable variables L, h, 1, H, and q, serves as a basis for 7 gg 
most of the work to follow inthis section andis used in its exact form, as Eq. x | 
in derivations. However, when it is to be used for computing Lin terms of — 


known values of H, 1, h, and q, it is much — convenient towork with a mod- : 


ified form that can be as follows: 


Eq. was derived consistently dropping all powers of r the third. 
It is very easy to use and is sufficiently precise for most guys. In fact, it is 
identical in form with a formula presented for the same purpose by | c. , Mor- _ 
: Tis and S. T. Carpenter,3 however, since the authors of that work were using 
properties of the parabola instead the their formula expressed in 
the symbols « of this study, 


is seen to be in error by a factor, cosy for the amount added to the chord | 


_ Although Eq. 16c is adequate for normal computations with guys, it may not | 


be precise enough for L L for c: cables with excessive For such 
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cables” it may be necessary tore retain n additional terms in the expansion of sinh r. 
C7 retaining powers ofr including the fifth, one obtains: _ 


(16e 


2 
i Effects « of Changes in Cable Parameters. —Each time a tower is loaded or 


 guyed tower, one must. have available formulae that ‘accurately relate | these 
. ‘ parameters and their changes. The three main classes of these changes are 


(1 li changes due to a load, affecting H, H, L, and h. 


2) h and L change due to a variation in temperature, ‘affecting » 7 


‘relation 


The partial in Eq. are Eq. 9 and sim- 
3. 


Note that inderiving Eqs. 18, 19, and 20, hyperbolic functions were Sita ne 


_ can be obtained by causing one to partially cancel the other. For example, if 
four terms of the series were instead of three, Eq. 18 would be written: 


— 
-« 
undergoes a change of dimensions, some of the Cable parameters must vary a 
 . 
— 
— 
_ — 
essary to computate quantitative values for the three maj 
above. Cases 1 and 2 may be derived from the differential 
ar 
= 
was frequently used. If care is excerised with these two approximations, 
"errors willnot be accumulative and in many instances an increase inaccuracy 
— 


The final form presented as as 18 can be r 


= aH _ 12 HS 


Me 
Here, cancellation of the involved in the two 


where the variables with the subscript the the initial values before 
change, or AH =H - Hp; Al=1- 1p; Ah = and AL=L-Ly 
Since r contains H, it is not constant during the integration of the preceding _ 
expressions, however, according tothe mean value theorem for integrals, there © 
«exists an intermediate value of r, denoted rj, Such | that r can be replaced by ay 
ry and treatedas a constant during integration. The values of the are 
not sensitive to small changes in r as r2 <1, If the initial and final values of _ 
r are known, rj should be made the arithmetical mean of the two values, but | 
when the initial and final values are not both known, it is necessary to a i 


p ification, yields, 


we 
tine 6 
oe K<1, Eq. 25 is easier to use in the form 


= 


_ The terms AL and Ah are seldom known. Even in case 2, AL and Ah are 
not completely known, because a change in H causes a change in L and h dueto 
% fact that the cable and tower are elastic. Therefore, it is necessary to be _ 
able to express AL and Ah as functions of AH in Eq. 26a before it can be used. — 


> 
> 


— 
— 
our 
— 
— 
— 
BUESS De Made, te COMpulalion Call Ue Pepedated usilig the TT's 
indicated after some algebraic sim- 
| 
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(T - To) 1 


e 


which Ag Ac and Ec are the area modulus, o of respectively, « of th 


cable, and T and To are, respectively, the finaland initial cable tensions. Note 
 thatasr 0, the result of Eq. 27 approaches the result obtained by consider-_ 
the cable to be a straight two-force member. 
_ An expression for Ah could include factors other than temperature; namely, 
_compression of leg of tower due to change in cable stress or decrease in length 
_ due to bending. - These factors could be considered, but since the effect of Ah 7 
_ usually smallas compared with other effects, the computations are not war- _ 
ranted. _ Therefore, it isthe writer’s recommendation that the cable stress and 
‘ee factors be approximated by assuming that the change in cable stress _ 
compresses or permits elongation of amember which has a cross-section equal ae 
_ the area of the leg at the level of the cable connection and a seninn equal tc to * 


‘From Eq. 11, the vertical component of ¢ of cable force, V, is 


Some: r+ L) 
esch l-<sr 2) an 


obtained by treating A cable as ‘atwo- force member, but the approach is. from : 
_ below instead of from above as one may expect, paso 
cs ‘Now that AV is known | in terms of AH, one can make use of the r relation . Ah = = a 


Ay and 1 E} are the area andelastic modulus of the tower leg at the level 


which can be with Eq. 27 to an expression for K that will account 


‘a cable stretch and compression of the tower leg. _ For example, AH can be : 


— 


‘determined for a given value of Al using Eq. 26a, in which is — as 


8 sti 
| 
— 
| 
— 
— 
As will be subsequently shown, this is easy to work with numerically. 


ee the formulae developed thus far to the computation of the total resistance | : 
offered, by allguys attached tothe tower at one level, to motion in any direction. a 
“Inspection ¢ of the formulae presented for the solution of a single guy, shows that _ 
the relation of guy resistance to tower motion is not linear. To evaluate the 
_ resistance to a given motion, one must employ a ‘trial and error’ procedure | 
for solving the equations. ‘This ‘ ‘trial and error’ approach is facilitated by 
having formulae expressing Alasa Senetion of AH for each cable instead 


_ By substituting the equations that express AL and Ah in terms of AH into’ 


35a, the following working equation is obtained: = 


Al= 


4 ‘miei of the us use of Eq. 35b wills subsequently be | presented to de- 
= the tower movement necessary for a group of guys to develop a given 
3 total resistance in a specified direction, 


ve 


= 


met In a detailed study of the efficiency of various cable arrangements, it is nec- ai 
essary to investigate the effect on the tower of wind from different directions. i 
han investigation is complicated by the fact that the relation between 


— 
— 
— 
nd 34, Eq. 32 and integrating, one obtains 
14057) 
‘ 
— 
— 
| 
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oe held constant in the derivation of Eq. 27, the mathematical model is elas- 
# €6=>_ tic but nonlinear, that is, loading and unloading occur along the same curve and 


+ 


A tower guyed by straight two-force ‘members has a mathematical model 
‘that is both linear and elastic. The direction of tower movement, in this case, | 
- depends only on the direction of the wind force and is not effected by its mag- 7 


nitude. _If the force is applied ineither of the orthogonal * principal” directions, 


the movement is parallel to the applied force. If, in addition, the “spring con- 
in these two directions are equal, as inthe case of a guying system with 


two axes of symmetry, the tower movement is parallel to the applied force for - 


In the case of a tower guyed with cables having nonlinear sag effects, the — 


at ‘direction of the tower movement is a function of the magnitude as well as the : 
. direction of the applied force and the movement is along a curved path as the 

magnitude of the applied force increases. The tower movement is straight and 
the applied force only if the force is along axis of ‘symmetry. 


: ‘Fig. 2 shows three common guy arrangements employed in the construction 
i tall towers. The loading is at an angle, 6. In none of these examples would = 
the 3 movement of the center of the tower be along the line of action of Rfor a 
general angle, 8. However, in any of the three cases shown, if @ equals 0° or 
60° and cable properties are symmetrical, a condition of total symmetry exists. 
and the tower movement will be along the line of actionR. eine 
_ Before proceeding with numerical examples, it is necessary to discuss a a 
"method for determining Al for a guy whose horizontal projection is not parallel 
to the tower movement. If the cable projection and line of movement intersect 
at some angle, ¢, then for atower movement, y, that is very smallas compared 
_ with the horizontal cable projection, 1, one can usually write Al = y cos ¢, that 
is, Al is taken as the component of y along the cable project however, the — 


component of y perpendicular to cable projection, y sin 9, does" cause 
increase in the length of the sin 9)? and for 


cise computations, it may be necess ry to c ae 


‘ 
i 
| 
7 
— 
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Numerical Examples. — Assuming ‘the constants: s: 1 = 950 he 
1349 ft; Ac = = 2.42 sq in.; Eg = 24. - ksi; a= hy 0 pif; A = 26.4 sq in.; E} “4 
30,103 ksi; Ho = 85.3 k; ro = 0.1211; = 0.0147; and L= “12+ he = 


1649.94 ft; then 


1649.94 + 1.34 + 0.003 - 0.0005 = 1651.28 ft ‘ 


> developed in ‘previous peal . The . physical constants ¢ given are 4 
_ those of the top cables that support an existing tower. = 

P ppo g 


. 1) Find the effect on Hy of a temperature drop of 40° F. Eq. 26b can be 2 
used for the trial solution to this problem, where Eq. 25 can be applied. The 

_ change in H will be small as the decrease in L and h are | cancelling « effects, sO 


will be taken equal to r?, that. is, rs = 0. our. Taking the co coefficient of ex- 


pansion for st steel, a, equal to 6.7 ?. 10° o/ee, and using the Eqs. 27 and 30, one 


40°) + +; 2) AH = -0. 442 + 0.0489 


- H = -0.362 - 0.0048 AH 
ae) the expression for Ah, t the factor 2 was inserted to account for the fact that 
“two cables are attached to each leg, as in Fig. — 
‘With a = 0, the working formula for K is 
K= 4. 53.10- = 0.109 [1- 0.361 AH] 


Try 4H = = 1.0k, the K = 0697 


sa) (2 (1.0553) = 1. 0k 


y AH = 1.20k, then K = 0.0616 and AH = 1.14k. aaa 
ow iment to make the difference between the trial value and the re- 
zero, one obtains AH = 1.167 which yields AH = 1. 167k substituted in- 


oo 2) Find the additional ene of an increase of q q due to 21 in. ice on the er entire | 
a e.. _ Here the increase in q will be dealt with separately and then the result will | 
a ‘be corrected to account f for cable stretchand tower compression. If L, h and] 


— 
= 
iim 
— Ly E+ — 
— 
4 
~ 
ig 
— 
= 
he equation being solved for AH also contains AH. so that a trialso-_— 
 &§ 
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Ah = -0.0176 - 0.0048 AH; or K = -0.0422 [AH + 0. a: 
s a net AH = 3.5k, then K = -0.149 and oa 
| 0.112 + 0.01 0.002 
or net AH = 6.3 - 2. 8,= 3.5k as ‘assumed. 
‘Thus, with ice on the cable, H = 40 kips and r= 0. 125. Then the vertic veutien! 


[neo coth r + = 


3) Find the effect of Al- = 2. 5 ft. 


. Here the change in H is probably going to be | on » the order of magnitude of HD 

go that Eq. 26a should be used. _ Since all the changes inh and L are functions ; 

of 4H aaa the definition for K given in Eq. 31 should be used. iw dot 


0.4 al - 0.0892 ¢- -0.713 an] 


Guess AH = = 15. 1, 0. 011 and 


ve: 


4) *F Find the motion of ‘the necessary for the guys to a siren 
™ re g the forces 

and moments to the tower motion. . This example will be based o on n the same in- 
conditions as used in example 

» Fig. 3, all six cables have identical properties, but are numbered as lor 

(2 to distinguish between the cables whose horizontal projections are parallel | 

to the line of action of the force and those whose horizontal projections ter- 
sect it at 60°. The tented 1 and 2 inthe the following equations refer to cables 


Ad 
— 
must be decr einH. na Thi scom- . 
Bey “the tower due to the net chang resses due to its own ice we al — 
— 
— 
; 


~ 
/ 


cos 60°) (am) = 2(AH2 - 
1. 


sin n 60")? | 
The ‘recommended the cable stresses 


es caused 


4 by R, is as follows: (a) guess AH2 and compute the corresponding Ai, thet 


“is, AHy = = (b) ‘compute Al and Alg to check for consistency, that is 


s, 


Kips; and = = 0. 01562. Substi- 


+04 
then | 


uess = 14.5, then . = = 55. 
1= (0.01562) 0.024; 


om 
47 / 40 


+ (0.0891) (0. 9736) (- 44 ft 


— gg 
wel 
— 

= 
— 

| 


A eluding trial with, = 15. 3 nets a residual of (+0. 0.06). I Interpolating to 
make the residual zero, one obtains the correct results: =" = 15. 23 | kips; 
° 4H = = -14. 82 kips; and y = 5.00 ft for R = 60.1 kips. pales, 


As a result of the above tower movement, the vertical components of stress: 


in cables 2 increase and the vertical components of stress in cables 1 decrease. 


This causes siamo to be applied to the tower that may be computed as fol-— 
lows: 


M 


= 


or = = 85, = 21.53 and M = = 293.5 kip- 


Z Ih order to obtain information on the tangent to the curve whenet R and a 
s preceding computations are repeated for R = 62.0, giving AH2 = 15.91, AH 1 : 


= -15.09, M = 302.5 and y = 5.20 ft. Therefore, inthe range near R = 60.1 = 
the function R versus y for the guy system shown, f(y), can be written: f(y 
95 y+ 6. this the external moment may be writt written 


ia 


4AV9 + = 


DYNAMIC ANALYSIS 
Forced Vibration of Guy. Cables.—By y solving f the equations of 1 motion : for 
a vibrating guy cable, it is proposed to determine the variations in horizontal Z 
force exerted on the tower when the support point is executing harmonic mo- 
tion. For small vibrations, out-of-plane lateral vibrations have the same equa- 
_tions.as in-plane vibrations, however, for convenience of presentation, the work — 
will be devoted to in-plane lateral vibrations. 
_ The coordinates ¢€ and 7 in Fig. 4 are chosen along the. arc and iperpendicular — 
to the arc of the cable in its static equilibrium position. Summing forces on 


an element of cable in a direction perpendicular to its equilibrium position, 
that ts, is, in the direction n, 


one obtains 


— 
= 
— 


GUY CABL 
CABLES 


a where € and 7 arethe coordinates shown in Fig. 4, p isthe mass per unit length 
of cable, and T is the tension on oe ical —e- the cable is in its equilibrium — 


"position. For small cable vibrations, the second order coreee | on the left side 


the can be ignored, with the result 


, Substitution into Eq. yields, 


— 
— 


on n is assumed that all points on the cable are executing saanincundien so 


that the variables in Eq. 44 can be separated as follows: 7 = 7} ‘cos at where n a 


isa mode shape function of x only. Substituting into Eq. 44 
&. 
+ cot cosh (2rx +a 


By expanding the hyperbolic function into its series, it is wn to obtain 

a series solution tothe ordinary differential equation, Eq. . 45, however, the se- 
ries solution is qe difficult to work with and such a general solution is not — 
-. actually needed here. What is needed is the solution for a vibrating cable that 


is functioning as a guy, that is,a cable whose tension varies only slightly along 


its length. The problem of the vibration of a cable whose tension, or density, 
or combination of the two varies slowly along its length can be solved by the ¥ 
_W-K-B (Wentzel- -Kramers- Brillouin)” method that is used in wave mechanics — 
_ and other branches of physics. In Eq. 45, the variable coefficient appears in 


cable of varying density. © 


be: of the dependent variable, 7, so so that the form is that of of the equation | for 2 


“3 The governing differential equation for the vibration of a varying denaity 


46) 


not vary rapidly it in the range under ‘consideration, the tn according to 


Adapting this solution to Eq. , 45 and making use of the boundary condition, F 


[ cosh (2 rx + a,)| 


at 


which c is an integration constant to be evaluated by the t boundary condition — 


x = 1, or the solutionto Eq.44is 


nxt) 
“Mechanics Slater and N.H. Frank, ist Ed., McGraw-Hill Book Co., New 


.&§ 
| 
— 
g 
— 
i 


in which D(x) is fa function | of X as indicated by the integral in Eq. 48. This 7 
; = one not seem to be a standard function, however it can be evaluated ot 


wile 
(50a) ) 
a 


= Vcosh + “ary + (50) 
‘The boundary c condition at x=1 is 


7 in which y cos at is | the tower motion at the support point and tan Y1> is the 


slope of the guy at the tower when it is in its equilibrium position. ™ “sub- _ 
stitution of this condition into Eq. yields 


v1 [cosh (2r + ay 


Eq. 48a is the equation of motion for a guy cable that is being vibrated about 
its equilibrium position by a forced harmonic motion at its upper r end. While 
_ this result is of interest, it is not the objective sought. The problem is to de- 7 


the force on the as it the transverse y sin 


eos at, to the upper end of the: cable. In general, this force is equal to 72 


48a, 


(HA) sin 


cosh 2rx+aj;) ‘cot 
ay) /4 


cos (54a) 


4 
— D(x) = | ¥ cosh (2rx + a1) 
* Trey (58) ~ Tase)(24) * [256)(352) 
— 
— 


y ein y, | cot BDy 


4 


- mustexert on the cable to impart a lateral motion of y sin yj cos at. The tow- 


Py er also imparts al longitudinal cable motion of y cos yy cos at which can be 
_ accounted for by using the static equations developed inthe first section of this” 


The i at xX = the tow er 


4h y cos y, sin y 
By iateabeie use se of the Eqs. 18 and 20, one finds that the amplitude of the hori- 
force that the tower must exert on the the cable tc to 


cos 
Combining this om ‘the horizontal component of the previously stated trans- 
verse force, one obtains the expression for the additional horizontal force that | 
the tower must exert on the ; guy to impart a horizontal motion, y — a t, at its 


+2 (1 - .267 r tan y Ale cosat 


—_—, Example.— —The following will serve to illustrate the numerical 


use of Eq. 56 and its application to groups of guys for finding the total dynamic 
‘Tesistance they offer to tower movement at a support point. ae. 


— 
—_— 
— 
a 
— 
>. 
| 
a 
56) 
— 


—_— same guy system will be used here as in the previous illustrative ex- 
amples. As shown in Fig. 5, the tower is assumed to be oscillating so as to 
cause symmetrical changes in the guy stresses. The line of action of the re- 
sultant forces is parallel to the horizontal projection of the cables denoted . 
fa _and intersects the projection of the cables denoted 2 at 60°, that is, the 1 cables | 
‘ are experiencing only an in-plane forced horizontal vibration of y cos at and a 
the 2 cables are experiencing both | an in- ee motion y cos 60° cos at and an 


— 


ih Eq. 56, the first term accounts for the transverse component of cable 
motion and the second term gives the static force necessary to impart the lon- : 
gitudinal component of motion. Because only the first term is used in express- 


ing the force necessary to impart out- -of-plane motion, it is convenient to re-_ 
oy write Eq. 56 using compacted 1 notation for the two terms, that is, 


‘FIG, 5.—SYMMETRICAL VIBRATION ¢ OF ' ‘TOWER AT SUPPORT POINT 


2,0) 3 


Bese ‘it may be seen that the resultant force inthe direction of the line of move-— 


of the tower is =" 


) 


r is vibrating abo oe a vertical equilibrium position, then T(r, 


— 
= 
an 
ay 
> 
and Ly ,) tie same rer ait 


as 

) 

+3 L(y yeosat .... 

ki similar expression may also be obtained for the dynamic moment caused 
by the changes inthe vertical components of the guy forces as the « support point o & 
oscillates. Only in-plane vibrations will cause changes in the vertical tee 
ponents of the guy f forces, use of Eq. 29, one 


a example, t take = 40. 0 = +1255, h= = = 1349 ft; 1: = = 950 
yy = 1 = 1.0255; and q = .01055 kips per ft. Then 


02242; = 1- 


1.3186; a) = 1161 a cot ( 


a 


horizontal of cable force 


ay 


4 
= in == (1 -= 12) )) Yoosat .......(59>) 
The moment, %(y), may be expressed as 
— 
= 
= 
Zz cosh a, -1 = 0.75497; z 
259; 
— 
— 
= 53] y cos at 
— 
— 
id 


CABLES 


1 
cable parameter equal: to: 
arameter 


H 


areas and moduli of elasticity for cable and tower leg, all a 


“normal and arc coordinates 


mass of cable 


special 1 integral function for dynamic solution 
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_ LATERAL IN INELASTIC BUCKLING OF TIED ssc 


Chin Fung Kae" 


- elastic limit is exceeded. The behavior of the arch in the light of existing in- 

- elastic buckling theories is discussed and a theoretical method for computing _ 
the inelastic buckling load is produced. 
The theoretical analysis is followed by anexperimental study with two main 

_ objects - the first to verify the results of computation, and the second to examine 

- the validity of a method of plotting suggested by R. V. Southwellfor certain types” 

of elastic buckling, when extended beyond the original conception to problems — 


» of inelastic 


The problem of buckling is one of stability and not of flexural stresses, and 
. a phenomenon peculiar to perfect structural members that are subject to- 
compressive loads. The general basis “of structural design is the evaluation 
of stresses onthe tacit assumption that equilibrium exists between the external 
_ and t the internal forces. The permissible stresses, therefore, determine the 


Note, —Discussion o open until June 1, 1961. To extend the closing date one month, a 

written request must be filed with the Executive Secretary, ASCE, This paper is part — 
of the copyrighted Journal of the Structural Division, Proceedings of the American So- 

of Civil Engineers, Vol. 87, No.ST 1, January,1961, 

I Senior Lecturer in Engrg., Univ. of Malaya, Kuala Lumpur, Malaya, ee 

“The Lateral Buckling of Tied Arches, by Ww. Godden, Part Ill, Proceedings, 

~ “The Design of the Uniraced Stabbogen Arch, ,” by Chin Fung Kee, The Structural © 


Engin ineer, Journal I, Struct. E. 


— January, 1961 sti 
_ 
A 
Fer 
— 
a 2 _ This paper is an extension of previous works.2,3 It concerns the same struc- a 

— 


January, , 1961 
and of the members. The strength of certain mem- 
bers, like the ideal column and the ideal stabbogen arch, however, is not gov- 
i: erned by the avoidance of a certainworking stress or of a fraction thereof, _ 
by the critical load at whicha peculiar | state of unstable equilibrium occurs. 7 
This critical load is called the buckling load. The inception of buckling is the Si 


i. result of a sudden breakdown of the internal resistance of the compressive mem- - 


€ ber regardless of whether or not the maximum stress exceeds that of the ma- a 
terial at the proportional limit at the instant of failure. . if it takes place ° when 

i" the maximum stress is below the limit towhich Hooke’s law remains yoyo 

the buckling is said to be elastic; if it occurs when the proportional limit has 


- _ been exceeded, the phenomenon is referred to as inelastic buckling. The latter 
< will arise in the case of “stiff” or non-slender arch ribs and in arch “ibe 


of certain metals, like light alloys which do not have a camel modulus of 
elasticity throughout the working range of stress. 7 


_ Notation.—The letter symbols adopted for use in this papera are defined wher 
they first appear, in the illustrations or in the text, and are > arranged _ 


for convenience of reference, in the Appendix. 


THE STABBOGEN ARCH 


‘The shape of the arches in n elevation is parabolic and it is assumed that there | 
‘is no bending moment in the vertical plane XAY (Fig. 1) - an assumption which 
is only valid if the hangers are inextensible. arches are e loaded through the 


_ to take tension only and have small flexural rigidity. The 
: to which the hangers are > connected is assumed to be absolutely rigid. The load 


of the arch instead of the 


STABBOGEN ARCH 


q 
4 
— 
| 
— 
may be considered 2° by means ofa horizonial thrust ar 
— 
— end 86 
. ‘thrust at which the lateral deflected form of equilibrium of the originally ideal — 


will TIED ‘ARCHES 


arch, resulting from a small lateral disturbing force Un- 
like the elastic case, it is not sufficient to define it as the smallest load at which © 
_ the original straightform (in plan) of the ideal arch becomes unstable, for while - 
_ in the elastic case, the arch resumes its originally straight form with the re- 
- movalof H', no such resumption is possible once the elastic limit is exceeded. 4 
The distribution of stress over any cross-section ofthe arch rib will depend : 
on the history of loading of 4 If it is applied before any of the compressive — . 
_ stresses reach the proportional limit, the stresses due to bending and todirect 
_ compression increase simultaneously. On the other hand, if it is applied when — 
_ the compressive stresses have exceeded the elastic limit, the superimposition ~ 
_of the bending stresses will result ina reduction of stress over the convex side — 
(in plan) and an increased compressive stress on the concave side. Stress | re- “ 
versal or decrement occurs. 
Under the first process of loading, the tangent modulus theory applies. — a 
7 der this theory, it is assumed that the compressive stress, and therefore the 
tangent modulus is constant over the whole section under consideration. To 
= the second mode of loading, Von Karman’s double modulus the theory is Ss applicable. — 
For simplicity, the tangent modulus Et is assumed to be constant over the in- . 
elastic portion of the section. For the equilibrium of the internal forces, it can 
shown that E;/E is equal to the ratio of the cross-sectional areas to the left 
rd and right of the axis N-N. ad The axis N- N is the axis of bending and is the line a 
along which the stress has not changed | with the ‘Superimposition of the bending “ 4 : 
The reduced modulus is givenby 


+E Et Ip) 


in which and Ig are the moments of inertia ofthe cross-sectionalareas sepa-— 


rated by line N-N, with tothis lineas the axis, = 


Aen 


EVALUATION OF THE ELASTIC BI BUCKLING LOAD 


was produced by W. G. Godden.2 It makes use of the for 
bility, that is: that at the critical load, the work done by the external forces is — 
equal to the increase in the strain energy acquired by the the system. A deflected , 


form in is given by the equation 


in which Li is athe span of the arch. This is the eine expression 1 of the lateral - 
_ deformation satisfying all the boundary conditions of an arch encastre at the 
ends - - - also referred to = end conditions A in this paper. The ends of hye arch 


 oTFA is the shortening of the chord AB due to deformation of the arch inthe 
-XAZ plane, and H is the horizontal thrust, then the work done byH is \ H. ‘The — 

- 4 “Inelastic Column Theories and an Analysis of Experimental Observations,’ P by c. 
T. Journal of Aeronautical Science, May, 1947, 
of Metal Structures,” by F. Bleich, 1952, 


= 
4 
— 
| 
> j 
| 
— 
| 
| 
— 
— 
— 
— 
. 
— 


ncrease in in | the strain sien of | the s atlas as a result of the deformationis _ 


= 


into rele equal sections, Giving the 


~ tional to “a The same value of A can therefore be used for arches of gore 


pe ~ It is also not possible to integrate the right hand side of Eq. 3 directly, | and 
: it has to be computed in the same way as in the evaluation of A. Table lgives 


part of the calculation for an arch with & rise-span of 0.12. 
ai ‘THE VALUE OF THE ELASTIC BUCKLING LOAD ee 
8 ees the results obtained from Table 1 and Eq. 5 to Eq. 3, the value of 


_ the elastic buckling load for an arch with a rise-span ratio f = 0. 12, is given 


6936 — (0.3496 + 0.0057 y) L 


30.693 
~ (0.3496 + 0.0057 7) 


4 


In contrast to the inelastic buckling ofa column subject ‘ani toan axial com- 
_ pressive load, the s stresses and therefore, ‘the values of the tangent ane: 
= the ‘inelastic modulus of rigidity Ct vary along the length of the arch rib. 


= 


which B = EI an is the torsional rigidity. ereiore, a e critical load, 
— 
wee 
{ 
— 
METHOD OF CON LOAD BY 
cause E; and C; are not constant. 
CEE 


ue of the inelastic buckling voy then the average compressive 


is angle of inclination of the tangent 2 atP to the horizontal. “The 
value of fe is least aa the crown C of the arch, and increases to its maximum 
~ at the ends A and B. Hence, the value of the tangent modulus E; is greatest at 
C, and decreases to its minimum at A and B. The flexural rigidity E; L there- a 
fore varies along the arch. Hence, the problem of evaluating the inelastic 
; buckling load is reduced to one of computing the buckling load of an equivalent 
a arch with a non-uniform flexural rigidity. This equivalent arch cannot, — 
ever, be determined until the value of the inelastic buckling load h is known, 


since the flexural rigidity at each section 0 of the the ase rib is tamed on the the 


"TABLE 14,—PART OF COMPUTATION FOR 


x 
0146 


0.0428 | 0.0263 | 0.0018 
0 0465 | 0.0022 
0.0459 0067 | 0.0021 
0.0400 | 0. 0.0016 
0.0297 ‘1010 | 0.0009 
0.0151 0.1515 0. 0002 


M“ds = 0. 3496 


End condition A f= = 0. 12 
M denotes bending ‘moment in the arch rib about y's ‘sand 
© T denotes torque in arch rib about the axis 


of h, value of h can, however, be assumed first in order to determine 
= the value of E; at each section of the equivalent arch, the buckling load of which > 
— - is then computed 1 by the Godden method. it the value of f the buckling load so ob- 
a _ tained is equal to that of h, then the value of h is the r required buckling load. — 
_ In assuming a value of h, one is guided by the limits: (1) h is less than H,;,, 
7 since at no section of the archis the value of the tangent modulus greater than © 
_E; and (2) h is greater than the » buckling load of an arch with aconstant flexural 
rigidity Ey 1, if Et is the tangent modulus corresponding to the average com- 
pressive stress (Her cos + V sin 0)/A at the ends AandB. 
The method of computing the inelastic buckling load by successive approxi- — 
mation may be outlined as follows; 
Pa... Step I - The Upper Bound.—An estimate of the inelastic buckling load is first — 
obtained by substituting the value of the elastic modulus E of the material in 


— | 
— 
— | ds/ax | M*ds/dx | 
= Section | | | — 
2 x0 4 H“ x0 0.0000 
— 0 | Hx0 | H 1,0929 >< 
q 0.2100 | 0.0074 | 0.0441 | 0.0001 0.0765 | 
a 2 1.0498 | 0.0824 .0007 
0.2410 1.0284 | 0. 22 
‘| 0.1623 | 0.0024 | 0,00 
| -0.0821 —1,0072 | 0.0446 | 0.0016 
} 1.0000} 0.1716 | — 
4 
— _ 
4 
| 


in which K is a on rise-span ratio the arch.? Eq. 7 
‘ gives the upper bound for the approximation of h, since the value of the e tangent 
modulus at every section along the archis less thanE. __ 
_ Step IT - The Lower Bound.—From the stress-strain curve, the tangent mod- i 
ulus E; corresponding to the compressive stress (Hoy cos @ + V sin 6)/A is 
found. (A is the cross-sectional area of the arch rib.) This value of E; is lower 
than the tangent modulus at any section of the arch, and, hence, the lower limit | : 
of h is obtained by substituting this value of Et for E in Eq. 7. This load -. ; 


be denoted by The inelastic should a value lying between “7 


0,00000H? | 0,20729H 1,10869H 
|0,01457 0.00000 0.17758 | 1, 
(|0.04750 |0,00005 | 0.14856 | 1.07706 
0 | | 1.06283 
0.00069 | 0.09754 | 1.04984 
0.00188 0. 05604 | 1.02844 
0.00220 03920 | 1.01980 
0.02531 | 
0.01430 | 1.00720 


100822 
1.00080 


1 1.00000 


‘™ his s table s applies to both the e double an andthe tangent Ss 
(The computations" indicate that by substituting the tangent modulus corre- 
sponding to the compressive stress Hcy/A instead of that corresponding 
cos 6 + V sin @)/A, for E in Eq. 7, a closer approximation is obtained.) 
Step III,—Using the values of Hi, the values of the tangent moduliat the vs — 


- ous cross-sections along the arch _e be found. These values are then used in ? ps 
the summation Gx a closer approximation on H2 


is shown in the accompanying is hardly re- 

pent Step II more than twice before a close of the inelastic 
os THE INELASTIC BUCKLING LOAD BY THE TANGENT MODULUS 
_ In the method previously outlined, since the tangent ‘modulus is used in the _ 


‘summation - he inelastic load so obtained is 


— 
(0.4329 | 0.9014 | 0.4800 
| 0.9150 | 0.4400 
| 0.9285 | 0.4000 
| 0.3387 | 0.9409 | 0.3600 
0.3048 | 0.9523 | 0.3200 
02695 | 0.9629 | 0.2800 
0.1960 | 0.9806 | 0.2000 
4 | 0.1582 | 0,9874 | 0.1600 
01192 | 0,9929 | 0.1200 
«(0.0802 | 0.9968 | 0.0800 0.10132 | 0.00088 0.00642 
| 0.0401 | 0.9992 | 0.0400 0.15160 0.00023 | 0.00160 
6 | | 10000} | #£=| o | 
| & 
— 
— 


the modulus load | Er i is instead of the | double modulus 
‘load hy is obtained. The former is obviously the quicker toevaluate as the tan- 
gent modulus E; can be directly read from the stress-strain relationship ——- 
‘the any given stress, while the evaluation of Ey is more involved. The inelastic 
buckling load obtained by the double modulus | theory is larger than that eek 
by the modulus theory as E, is always grenter than Et. 


— THE VALUE OF THE INELASTIC — OF RIGIDITY CE 


p between the | 
tangent modulus of the e tangent modulus Et. Different values 


(a) Ce =C C - This” expression is used in the theory of local buckling 
of plate elements. The term C is the elastic modulus of rigidity and Ey the 
(b) C/E - This relationship the treatment of 


(ce) Cp = =ErC/E - This is used by W. Ramberg and S. Levy.& Levy.® ‘The The terms 
Cr and E,. are the e uivalent moduli under the | Double Modulus theory; and re 
(d) {/E+ rather complicated expression is used 


it was decided to apply the relationships given by (b) and (c) in computing © 

the inelastic buckling loads under the tangent modulus and the double modulus _ 

theories respectively, since in each case, it would make the ratio of the lateral 

flexural to the torsional rigidities yt constant throughout both the elastic and 

s inelastic ranges. A comparison of the values of Yt derived by (a) and those 

tt either (b) or ( (c) for the model arches eC in the experiments, is as follows: pA 

( (b) or (c) | By (a) 


«1816 
316 1175 


The expressions M: and T are independent - of y, and the effect of y on the © 
buckling load comes only from the evaluation of the right hand side of Eq. 3. 
Lo 4 Since the value of © T2 6s is smallas compared to that of 2 M2 5s, the effect 
- of small variations in the value of y on the buckling load is silane is 


“On the Strength of Columns that Fail by Twisting, » by 


Science, Vol. 4, 1937, 


— P 
+ 
ve 
2 
— 
x: 
“Instahility of Extrusions Under Compressive Loads by W. Ramberg and S., Levy, 
. Lundquist, Journal of A 


shown by the results for a a rise- “span ratio of 0. 


316” 


1100 
| ‘THE V VALUE OF THE VERTICAL Sk SHEAR 


on 


rk: =a 


because there is no no resulting bending moment in the KAY But 


7 In the ‘following examples, the inelastic weatinas: loads are competed for 
model arches for which the span is 40 in. The tubular section has an external | yi 
7 diameter of 3/8 in. and a wall thickness 0.022 in., the cross-sectional area, A, — 
Se is 0.02442 ‘sq in., ly equals 0.0003814 in.,4 the value of E is 11.24 x a a 
and yt =y = 1. 316 . The values of Eyare read directly from Fig. 2, which gives 
the between the stress and the tangent modulus for the material of 
which the model arches are made. The value of n is obtainedfrom Fig. 3, which 
b _ gives the relationship between n and the ratio Et/E for a tubular section. . The 


with n for a tubular section of mean radius r and wall thickness t. _ eee pa 
‘The integrals are computed by the mid-ordinate rule by dividing the span in- - 7 
to twenty- -four ec equal sections. The computations were done ona calculating | 
No.1 - - The Inelastic Buckling Load the Tangent Modulus The- 


.—The rise- -span ratio f was 0.12 


. ithe of Ij and Ig are taken from Fig. 4, which gives the variation of 1; and 12 


— 
= 
in which f is the rise span rati 
— ii 
_ The vertical shear force at P is 
average compressive stress at P is (H cos V sin ; 
a 


/A = = 231. 02442 - = 9459 
ia From re. 2, ‘the value of Et for this stress is 9. 45 x 108 pe psi. Tevetore, 


Hy = = (86.1 x 9.45 x 106 x 0. 0003814) /402 = 194, 01 a ae 


Hy = 194, 194.0/0.02442 = = 7,942 psi 


Hy 202. 0 » (from ‘Table > 3) 


Hg 199. 6 Table e 3) 


“hy = = + (He + Hg) = 200.8 Ib 
Example No. 2 - The. Inelastic’ Buckling Load by the Double Modulus The- 7 


ory. —The rise- Span I ratio is 0.12. the final is (Table 


211. (trom Table 4); 


es 


hy, = 212.45 lb 
_ Introduction. —The objects of the experimental study were 2 to verify the com- 
puted values of the buckling loads, and to investigate the 7 


Though Southwell postulated his method for the case of the pin- 1-ended strut,8 
_ the method was found to be valid in predicting the elastic buckling load of the 
-_ _ stabbogen ; arch. ‘Southwell indicated thathis method would only apply within t the 
- elastic range, but C. T. Wang? has claimed its validity for predicting the in- 
elastic buckling loads of columns. Wang applied the method to the experimental - 
_ data obtained by Von Karman, Gerard and those of Horsfall and Sandorff, and it 
was found that the results gave good linear relationships | between 6/P and 5, 
in which 6 was the deflection resulting from a compressive load P. 
_ deemed that Southwell’s method might have wider applications than atfirst re-_ 
alized, and that it be the prediction of buckling loads instruc- 


in elastic buckling was employed. 2 The. only ‘Modification was i in nthe fixation of 
the en ends of the model arches to the end- plates by means of tappering plugs. 
_ The Model Arches.—™Model arches of commercially pure aluminum were 
used because this material 1 has (1) a comparatively - large inelastic range, and 
(2) a low e elastic limit as ‘compared to that of structural steel, and its s employ- 7 
“ment entailed very few modifications to the existing apparatus. = ~~ 
_ The section was of 24 gage with an external diameter of 3/8 in. The span of 7 
_the arches had to be reduced from 60 in. (the span 1 of the steel arches used in 


the experiments on elastic buckling), to 40 in. in order that buckling occurred 


beyond the elastic limit. 


8 “On the Analysis of Experimental Observations In Problems of Elastic Stability,” # 


by R. V. Southwell, Proceedings, Royal Soc., London, Series A, , Vol. (135, 1932, ‘ond 


— 
4 — 
Step 
— 
q 
. 
4 
= 


The model were parabolic in n elevation. Seven wire-hangers, spaced 
at equal intervals, were used. Each wire was wrapped round the boom of the 
_archand the ends were twined together | before clamping them to the trolley. > 
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FIG, 3.- —THE VARIATION OF n WITH E aie 


: in the tests was to observe lateral deflections with the corresponding end thrust © 
; va measured on both the proving rin rings. To ascertain that the experimental obser- 
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vations ma made corresponded t to stresses ‘beyond the. elastic Limit,, a quick stress 
a analysis and observations toensure the presence of permanent set were made. an 
ae The Southwell plot was obtained by using the values of (z - - Zz) and (z - - Zy/ 
i, (H- Hj) instead of those of z and z/H. _ The term (z - 2, was the amount by 
: = the lateral deflections were increased \ when the load was increased from 7 


= 
0-03r 


Al. 


: ane Southwell Plot. —The Southwell plot is is fairly linear if f plotted overa range 
"oof stress within which there is a negligible difference in the value of Eee Fig. <a 
the results of tests No. (18 and No. 22. For those plotted over a range 


cates the registration of a measurable ee in the value of Et. — 


— 
— stresses beyond the elastic range. 
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The Value Of The Inelastic Buckling Loads.- —The inelastic buckling loads as as 
assessed from the Southwell plots are only approximate and may be sufficiently - 
accurate for design purposes. Table 5 gives some of the experimental results 
in comparison with those evaluated by calculation under the tangent modulus | 
and the double modulus theories. The percentage error given in Table 5 is that 

based on the buckling loads computed on the tangent modulus theory. = 7 

As is indicatedin Fig. 7, which gives | the plot of the results, the experimen- 
tal results for the buckling loads are more closely related to those obtained by : 


“TABLE 5. —COMPARISON ¢ OF RESULTS 


Double Modulus Tangent Modulus 


Experimental Percentage par 


7 


“= 


209.6 
5 


7 
| 


g 


° 


FIG. 7. EXPERIMENTAL RESULTS AND THEORETICAL CURVES 
the tangent modulus es those evaluated — the double modulus theory. 
This is obvious because as a result of unavoidable eccentricities such as im- 


perfect alinement of the arch rib and a non-axial application of the horizontal 
end thrust, an arch under test will deflect laterally as soonas a load is applied. 
This results ina simultaneous increase in the stresses due to both bending and 
to direct compression complying thereby with the ential of the tangent 


modulus theory. 


| 
= 

oS 
— 
" 
if 
— 
ta 

3.71 
— 0.00 


nr _The method of successive approximation gives the value of th the 
"inelastic buckling load accurate e enough for design purposes. | The inelastic 


by the double modulus po dhesael and is, therefore, the safer load to use in 1 design. - 
_ The tangent modulus load, apart from being the easier one to evaluate, , gives a 


“i 


_ Southwell’s experimental method is only useful for predicting an ap approxi-— 
_ mation of the inelastic buckling load if the arches are of material which has a 
comparatively small variation in the tangent modulus. Because of its larger 
scatter, the Southwell plot lacks the advantage of serving as a check against — 
possible experimental errors due, for example, to slip inthe wire-hangers ot 
to local instability. The test specimens have to be of greater perfection than 
is necessary with those that buckle within the elastic range in order that the | 


» Soe load reached is a large proportion of the inelastic buckling load. 7 
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APPENDIX. NOTATION 


sectional a area of rib; 
Flexural Tigidity; e Bo = Ely; 


‘Torsional rigidity; 
Modulus of of 
Reduced modulus; 


— 
Tangent modulus; 


= = for the elastic range; hence fo for the inelastic range; 


zontal component of thrust in the arch 
“Critical v value of ‘Hin the elastic on that is the elastic 


The inelastic horizontal buckling load; 


i. ‘The inelastic horizontal buckling load by the double modu- 


The inelastic horizontal Gucieing load by the tangent modu- 
theory 


&§ 
4 L 

q 
7 
— 
— 


" _ Bending moment in the arch rib about the the axis 3 y's 


Displacementof the axis N- axis Y-Y 
(double modulus theory); 


“Mean radius of | > 


Length « of the arch rib center line; 


origin A measured the center 
line ofthe arch; 


7 Torque in the arch rib b about ' the axis x’; 


= Thickness s of tube; 
= = tan” 


vorttcs’ force on on chord AB; 


Rectan lar co- -ordinates; 


om Axis tangential and normal respectively to mnie line 
_ of the arch in the plane Sa; and 


ular "to the XAY 
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— 
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SYNOPSIS 


A method of analysis of elastic structures tant avoids: the solution of large 


_ Among the many ny refinements proposed for the | analysis of statically inde-— 
terminate structures, several methods have in the past been suggested that 
the solution of large systems of simultaneous equations. 
Among early investigators, Mueller-Breslau (1897),3 Mueller (1907), 4 and 
ae (1921)5 proposed the application of redundants in such combinations — 


as toresult ina diagonalized set of compatibility conditions, in which each equa- a ; 


tion contains one unknown explicitly. This method is generally applicable 

_ Note.—Discussion open until June 1 1, 1961, To extend the closing date one month a 
_ written request must be filed with the Executive Secretary, ASCE, 


This paper is part — 
of the copyrighted Journal of the Structural Division, arate of the American weak 
_ ciety of Civil Engineers, Vol. 87, No. ST 1, January, 1961, 

4 Research Asst., 


Brown Univ., Providence, R, 
2 assoc, Prof, of Civ. Engrg., Univ. of Colorado, Boulder, Colo, — 


4 “Zentralblatt der Bauverwaltung,” by H. Mueller-Breslau, 1897, p. 513, 


“Zentralblatt der Bauverwaltung,” by H. Mueller-Breslau, 1907, p. 23. 
5 “Der Eisenbau,” 1921, p. 305, 


— 
— 
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AS 
‘structures of any type or degree of indeterminacy and, in special — 
duces to the well known elastic center and column analogy 
&§ The underlying theory is presented, its use is indicated by examples, 
physical meaning is discussed, and its relat 
— 
— 
en 
i 
_ 
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and simple ‘set of operations suitable for | tabular or computer 
7 In cases in which either the degree of indeterminacy or the number of loading = 
conditions is large, this method a significant in the amount 
For the special case of flexural structures of the po circuit type, it can 
be shown that this approach is equivalent to the well known elastic center or 
column analogy methods of analysis. These more specialized procedures have 
a : received much attention in the literature of recent decades and are commonly 
— when many loading conditions are to be considered. __ 
_ It is the purpose of this paper to direct attention tothe more general formu- 
. lation of the method, to demonstrate its applicability to the analysis of various © 
«structures, and to establish with th the more e widely known pro- 


' OUTLINE ANI AND APPLICATION OF THE METHOD ; 
This « classical method of analyzing elastic structures consists of the reduc- a 
Hon of anindeterminate structure to a determinate form by removing anumber = 
of restraints (redundants) equal to the degree of indeterminacy and then reap- 


plying these redundants in such amounts as to satisfy the conditions of geometry. 


The set Eq -1 is solved for the X’s and the internal forces are then found by 

_ _The proposed method is an orthogonalization procedure, t by means of which | 
all interaction coefficients are eliminated, th thus leading toa system of equa- 


4 


— 
— 
“ite 
mu eous equations tor the redaundé 


TURES 
i To establish this method, “use e will be made of the principle of virtual 1 work fi 
because of its power and generality. However, the necessary equations are also 
a The principle of virtual work asserts that during a1 any y virtual displacement of 
- an elastic body the net work done by all the forces is zero. The internal and > 
external forces must be in equilibrium, and the virtual Sentnenet must be 
small and compatible with the conditions of constraint. 
4 For clarity of om, an application of the method toanalysis of a rather ; 


n generalized. ‘Consider the three times redundant, pin jointed truss of Fig. 
1 (a). The member forces S in this structure may be thought of as resulting | 7 Hl 
from superposition of the forces S, due to the external loads on the determinate § 
truss , and the panes le S,, due to “the redundants, where x, is the value of the ; 


redundant. We therefore write 


ié‘C 
— 
— 
— 8, +x, 8, 


s 


With supports, the deflections at points of supports var vanish, whereas” 
_ Consider next the truss of Fig. 1 (c) with a unit load as shown equilibrated — 
_ by bar forces sj and reactions that can be computed from. statics. Suppose this _ 
: force system in equilibrium is subject to a virtual displacement equal to the 
actual displacement of the truss of Fig. 1 (a) (that is assumed small inusual _ 
structural practice and is compatible with the ‘constraint. conditions). Then the 
net work done is zero. As load do work during 


The extends over members of the equation con- 
tains the unknowns xj, xg and x3. Two similar equations can be obtained by 
application of virtual work ¢ to the equilibrium systems of Figs. 1 (a) and 1 (e) 
subject to virtual displacements equal again to the actual ‘displacements of the 
truss of Fig. 1 (a). In this fashion, a set of equations of the form of Eq. 1 are 
_ set up, that can be solved simultaneously for the unknown values of the redun- 
ee It should - noticed that in this classical analysis the t truss systems of Figs. 7 
‘ 1 (c), 1 (d) and 1 (e) were chosen strictly for convenience. . In order to reduce — 
the resulting - set of equations to the form of Eq. es it is necessary to make a 
more appropriate choice of systems. A 
_ In particular, new systems of redundants are applied tothe determinate form 7 
of the truss that are combinations of the unit redundants of ‘Figs. 1 (c) tol(e). 
: These systems will be called combined unit redundants although i in generalonly 
one of the forces of the combination will have unit value. The bar forces due 
to the nth combined unit redundant will be designated 
fish. The first combined unit redundant will be chosene equal to the o one of Fig. 1 (c).. 


corresponding member forces are 
- 


member forces 


coefficients Cjj are constants to be determined. 


= 


Denoting by a the flexibility (compliance) of any bar (a= ‘A E) 
— 
i 
— 
| 
1 (a) such that the member forces are ray 1 
third combined unit redundant is a combination o : 
a 
oo 
if 


hing Eq. 3 superposition asserts thatthe member forces inthe 


= 


v 
x. s, +X, 
‘Here represents the value of the nth ‘combined unit redundant. The corres- 


ponding elongations are now expressed by 


S=a a(X, 8 +} 
Applying virtual work, 1 we subject each of the ‘unit redundant ‘sys- 
tems toa virtual displacement equal to the actual displacement of the real truss. 
As before, none of the external 1 forces 3 or reactions ¢ of the redundant systems = 7 
do any work during the displacement and we obtain 


s, 8 1S S.... (9 


e 


Th The summations again extend over all members of the truss. 7 
Ge Eqs. 9 constitute a set of the form of Eq. 1. —— to the _— scheme, a 


we must now make the cross terms vanish, that “7 


7 “This is a accomplished by y suitably choosing the coefficients Cij in Eq. 6. Sub- 


‘stituting the expressions from Eqs. 6a and 6b into Eq. 10a yields 


atistied if we choose ite 


7 
(9c) 
— 


The first term here is zero by Eq. 10a So we need only ll 


a $2 


> 2 
Beginning w with the known member forces Si. the constant Cor ‘can be com- 
_ puted from Eq. 11b and then forces Sg can be determined from Eq. 6b. Having 
obtained S2 the constants C3; and C39 can be compenen irom Eq. 11d and 11f. 

Bar Forces. 


— 
‘With these, the forces te can be determined from Eq. Sc. Deal combined u unit 


“redundants producing S;, Sg, and S3 are shown in Fig. ul 
are now reduced to the form 


“eee 


The multipliers Xp, are found a Eq. 12 and substituted into Eq. 7 to obtain 


the member forces. It iad that the Cmn depend only 


_ The second term, by » 
substitution into Eq. 10c gives 
Finally, substituti int give: 
. 
— 
yaa) 
| 


ANALYSIS OF STRUCTURES 
on the nature of the structure. For analysis of any other loading condition, only — 
‘Eqs.12and7arerequired, 
- Table 1 shows | the numerical c: calculation of the solution of the truss under > 
the load of Fig. 1. | For simplicity the flexibility a of the members has been | 
taken equalfor all members. 
_ Columns 1, 2, and 3 give member forces due to unit redundants as deter- _ 
‘mined by statics. Columns 4to 13 comprise the work of finding the coefficients © 
— Cmn and the member forces. due to the combined unit redundants. Columns 15 
_ to 17 enable determination of the X,'s, and columns 18 to 21 produce the final Ms 


Analyses for other loading conditions now st the work 


= of columns 14 to 21. 


GENERAL THEORY 


‘To formulate this method in a more comprehensive manner, , it is convenient 
- ® consider an m times statically indeterminate structure. Call the internal 
stresses s or S, and the corresponding strains 5. These stresses may be due | 
to axial force, shear force, bending, torsion, or any combinations thereof, pro-— 
viding all normal stresses at any point are unidirectional. 


The structure is made statically determinate by removing restraints. 


order that no external work will be done by interior redundants, removal of in- | 
. That is, in- 

ternal members are cut or hinged so that equal and opposite forces or moments _ 
can be applied at points infinitely close together. Therefore, the only external — 
work that can be done by a combined unit redundant system during its virtual 
displacement will be due to yielding of supports. 
a Unit values of the m removed redundants produce internal stresses $1, Sar 

. Sm, respectively in the determinate structure. From these redundants, 


a set of combined unit redundants that —" internal stresses Sj, So,.. 


Sm, is constructed according to the scheme; 4 


‘Beginning with Ca, the coefficients for e 
to 


— 
a 
— 
q 
se 
= 
4 
at 
— ach successive S, are computed ac- a 


q 


“bo 
ooo 


6 


| 
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Xm are by the principle 


| apie je ie je je je fe je’ ie je fe 
BC | 0,500] 0,250 |-0.707 | 0.250 | 0.125 ]-0,292 |-0.042 |-0.354| 0.030] 0.002|-0.147 |-0.003 
CD | 1,000] 0,500] 0 1,000 | 0.500|-0.583|-0.083/ 0 | 0 | 0.007|-0.294|-0.006 
DE | 0.500] 0.750| 0 0,250 | 0.375 |-0.292] 0.458 | 0 | 0,209|-0,147| 0,034 Ag 
FG |-0,500|-0.250/ 0 0,250 | 0,125] 0,292] 0.042 | 0 |0 | 0.002] 0.147] 0.003 
GH_ | -1,000 |-0,500 |-0.707 | 1.000 | 0.500] 0.583] 0.083 | 0.707 |-0.059| 0.007] 0.294] 0.006 — 
HJ | -0,500 |-0,.750 0.250 | 0.375| 0,292 |-0.458 0 | 0.209] 0.147\-0.034 
fo jo ljomlo Jo fo Jo Jo lo 
| -0,.500 |~0.250 |-0.707 | 0.250 0.292] 0,042 | 0.354 |-0.030 | 0.002] 0.147] 0.003 
CH | -0.500 |-0.250 |-0.707 | 0.250 0.292] 0,042} 0.354 |-0.030| 0.002] 0.147] 0.003 — 
DJ | 0.500 |-0.250 | 0 0.250 -0,292 |-0,542/ 0 | 0 0,293 |-0,147 |-0.041 — 
EK | 0,500] 0.750/ 0 0.250 -0.292| 0.458} 0 |0 | 0,209|-0,147] 0.034 
BF | 0.707] 0.354] 0 0.500 -0.412 |-0,058 0.003 |-0,.208 |-0.004 
€G | 0.707] 0.354] 1.00 | 0.500 | 0.gmo|-0.412 |-0.058 | 0.707 |-0.058 | 0.003 |-0,208 |-0.004 
0.354] 0 0.500 |-0.Mo| 0.412/ 0.766/0 | 0 | 0.587] 0.208] 0.058 
| -0,707 |-1.060 0.500 | 0.412/-0.648/0 | 0 | 0,420) 0.208/-0.049 
— —— 
which of proportionality be- . 
gives an orthogonal set of redundants such that — 
Call the in by applied 
loads So. E ‘ucture are 
— (16) 
— _ 
._ 
> 
+6,..07 
in which the 54 term denote: 
goon. The multipliers X1, — 


| 


0.643P 

0.150P 

0.028P 

0.750P} -0, J 0.112P 

.500P} -0. 0. -0.204P | - 
0.250P] -0.125P| -0.114P} 0.028P 


0.057P 
0.424P 
0.207P 
057P 


5 
we. 


0.139P |-0.250P 
0.139P 
0.136P] 0.047P] 0,250P 
-0.114P| 0.028P -0,250P 
0,061P| 0.224P| 0.3545 
0.279P | 
0.094P | - 
0.056P 


x 


= 


of virtual work to each the S1, subject 
- toa virtualdisplacement equalto the actual displacement of the real structure. 
In view of and the effect of support settlements, this gives 


A = 
2 


: The work resulting from yielding of the pth support by ar an | amount t Ap ii is included 


_ force or moment at the pth support of the : appropriate equilibrium system, and 
the corresponding settlement Ap- Calling the externalforce or momentat point — 
p for system Sp, rpp and for system Sp, Rpp it it can be seen that these are the — 


i) 
| 


(0.113 | 0,012 0,153P | 0.153P 

0.150 | 0.022 0.25 
0,557 0.310 |-0.250P] 0.125P 
— 0.188 | 0.036 | -0.250P] -0.125 0.284P | -0.216P 
0.012 | -0-250P] -0.125 0.096P | 0.096P 

-0.212| 0.045 | -1.060P| - o.os7P | 

— 0,159 | 0,026 | 0,354P| -0.250P 
— 1.763P| -0.081P |-0.081P 

— 

— 

— 

— 
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same for system Sy as for system 81, 01 or Rjp =Tip- For system S2, the Sup- 


port reactions are those for s2 plus C91 times those for $1, or Rap = Co} Rip 7 
system Sg they are those for system s3 plus times those for 
81 plus C32 times those for ‘So, or or Rgp C31 Rip + C32 R ap + 


allows for the successive determination of each Rn np” 
~ Each X is computed from Eq. 18, and the waieadl stresses of the structure 
are then obtained from Eq. 16. For other loading conditions applied to the same 
structure, the analysis requires only Eqs. 18 and 16. It should be noted that in 
_ Many analyses including influence line computations, the terms representing - 
‘support settlement and temperature : strains are customarily omitted because 
only the effects of live load are considered. = 
= Integration to obtain the internal work canbe simplified for structural mem-_ 
bers. Subject to the usual assumptions of structural theory, the state of —. 4 
-at any cross section of a member can be completely described by specifying — 
_ the bending moment, twisting moment, normal force and shear force that act 
at the section. . Call these M, T, N, and Vv, respectively. T The corresponding» 
elastic strains are also specified because of Hooke’s Law. 
5 Consider now the work dW done in a slice of infinitesimal thickness dl at a 
cross” section of a member in which the stresses of an equilibrium system Si ; 
are specified by Mi, M, The strains due to a virtual displacement 
system Sj ; and Vj . Integration over the cross section 


egration along the length of each lt inaania the entire structure 


PHYSICAL INTERPRETATION © 
The physical significance of the proposed procedure may be ‘Suggested bya 
consideration of the deflection curves of the three times indeterminate continu-— 
ous beam of Fig. 3. The statically determinate system selected is that of Fig. 
3 (b), and according to Eqs. 13, the combined unit redundants are those of Figs. 
_ Considering now, i in particular, the deflection curve of Fig. 3 (a), we see 
sj that the effect of selecting the coefficient C21 according to Eq. 14, is to annul . 
the deflectionat point 1. The deflection curve of Fig. 3 (e) due to the third com- a 
bined unit redundant is a superposition of curves proportional to those of Figs. 
3 (c) and 3 (d), and that due to the unit force at point 3, at which the amplitudes 
_ of the two first « curves, specified by C31. and C39, are such as to o annul deflec- aa 


tiosatpointstand2, 


= 
— 
— — 
an 
| 
= 
— 
— &§ 
; 
3 
— 
— 
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n 
~ 
oO 
~ 
22) 
< 
vA 


Wis 


: en of the deflected shapes due to the redundant states Si, So, S3,_ 
sk shows that, of these states, only state Sj will affect the deflection at point 1. _ 
The multiplier Xj is accordingly selected to annul the deflection at point 1 que 
to the load (state So). Of the other redundant states, only state So will = 


2 
the deflection at point 2, and XQ is — selected to annul the deflection 
at this point due to the states (Sp + X1 S 


_ The scheme is extended for addi- eg 
a tional redundants. Eqs. 18 provide for the ‘selection of appropriate values ot 7 
- Xp, and Eq. 16 for the superposition of the sates, tnt is also shown in the 
series of elastic curves of Fig. 3 ea ee 
it should be noted that the Sh determinate system So of Fig. 3 (b) ll 


which the work done by axial and shear forces can be neglected in comparison 
_ with that done by bending, and to show the relation of this method to the well 


7 known elastic center and column ‘analogy mé methods, consider a fixed ended arch 
structure, as shownin Fig.4. 


Make the structure determinate as shown in Fig. 4( (b). 
moment and horizontal and vertical forces 


T The unit redundants 


cut end, and the due to these are 


_ The combined unit redundants are ‘then selected ncneeding to me 013 | to give 1 mo- 


Cc 


The coefficients here are 14, as 


hae the analysis of continuous beams or frames might oftenbe to considerinternal 7m 
= RELATION TO OTHER METHODS 
— 


a ANALYSIS OF STRUCTURES 


fit- fet) /( 


The multipliers - are computed to Bq, 18, the final moments” 
are written by use of "Eq. . 16. It is, of course, not always necessary to co compute e 
the: integrals analytically. Sometimes, as for instance in the case of linear 
structures and moment diagrams, it is is more convenient to resort to graphical 
It is instructive to show the relation of this method to the well-known elastic _ 
center method applicable to this type of structure. In that method, the coupling - 7 
terms are eliminated by suitable choice of orientation and point of application 7 
of the redundants, namely along the axes with respect to which the quantities — 


as shown on Fig. 4 (c). The axes thus can identified as the centroidal 

principal axes xj, yj ofa conjugate structure the elements of which extend along 

the centerline of the structure a distance dil, and have a width normal to andl = : . 

centerline of 1/ElL The centroid of this “equivalent area is thus defined, 


terms of the axes x, y of of Fig. 4 (a), by tl the distances er 
& 


: sete the prnnaen of inertia of the conjugate structure with respect to the 
_ centroidal x2 and ‘V2 axes, and the denominator the moment of inertia with re- 


- 
& ‘> 

~ +& 
— 


t spect toa ‘centroidal y2 axis, So 


As 


dl 


This shows. that the elastic center method becomes a a special cas case of the pro- 
‘posed method for which all coefficients Ci vanish. ir 
_ _It also can be seen that if the x and y ‘axes through the point of redundancy a 
are taken parallel 1 to the principal axes + of the conjugate structure, then the co- 


_ particular, this occurs in symmetrical structures, thereby considerably 
a. For this symmetrical case, it will be shown that the computations for the 
oz proposed method, for the elastic center method, and for the method of column 
analogy become identical. 
a Here, Eq’ s. 23 reduce to 


=X +X 


— 
— 
— 
a 

— AR 


ANALYSIS OF STRUCTURES 


7 M, dl al aa x, 


2 Ei _J (x-x) “ET ET E (300) 


*) 


that will be recognized as the elastic center solution. Substituting th the usual 


‘symbols for the integral expressions, 6 this can be written as 


that is the column analogy form of the solution. vha 


‘CONCLUDING REMARKS | 


The present treatment isa specialization of the Schmidt Method of 


EA 


of stress ta structural analysis and elasticity has been by Prager. 
am In that paper, a geometrical approach t based on the representation of States of 


6 «statically Interdeterminate Structures,” by C. K. Wang, 1953. 

a 7 “Methods of Mathematical Physics,” by R. Courant and D, Hilbert, 1953, p. 4. 15 

- om, 8 “The Extremum Principles of the Mathematical Theory of Elasticity and Their Use 
in Stress Analysis,” by W. Prager, Univ. of Washington Engrg, Sta., Bulletin No. 119, 


4x” 
in which Mo denotes™the™ Mom due tie load onthe determinate 
ae — Here the coefficients Cy; and C3; are given by Eqs. 24a and 24b, andac- 
— 
a 
| 
| 
— 
| 
if 


any desired degree of accuracy. Morkovin and Kincaid have examined the sim- = 7 
_ plifications due to orthogonalization of stress states. 9 They, however, restricted 
their observations to symmetrical structures under bending moments only. 
fe For structures of a high degree of indeterminacy, some advantages are of- 
fe 


red by the Peach Method, 10 that uses matrix ‘methods to construct from the 


‘ling Eq. 15 has the — ge 


18. 
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Phil M. Ferguson! 


For concrete members, load- stress ci curves show irregular- 
— Gttes at cracking load and near ultimate load that the attainment of reasonably 
re predictions of member capacity — the use of ultimate strength 
_ Under ultimate strength theory, the designer must become familiar with load 


_ factors (and probability theory). Several load factor ‘proposals have been pre- 


~ over working loads. In addition, American study has always used separate load 
factors for dead and live load. 
The increasing use of strains andthe convenient use of the rectangular stress - 
block for concrete stresses seem to furnish the tools that the designer needs to 
develop a systematic “ theory” ultimate strength either her complicated or 


= sections. 


INTRODUCTION 


- willbe more efficient if he understands some of the background necessitat- 
7 ing the change to ultimate strength methods andis aware of — of the trends 


that must become part 0 of the ultimate strength picture. 


7 for the individual papers in this symposium. To extend the closing date one . month, a 
_ written request must be filed with the Executive Secretary, ASCE. This paper is part 
4 of the copyrighted Journal of the Structural Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. 1, , January, 1961, 
Prof. Civ. of Civ,  Engrg., Univ. of Texas, Austin 12, Tex, 
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a 

diversity of load-stress relationships why concepts based on working —— 
_ design are relatively inefficient tools. It indicatesthe probable increased role 

aa of the load factor concept in future design. Finally, it points out the necessity 
for a greater introduction of wamate Strain concepts into ultimate strength 


STRESS COM PLICATIONS | 
¥ ‘Ultimate. cial design is relatively new in codes, but it is not a new idea. 7 
Fifty years ago a form of ultimate strength design almost became the official © 
procedure for the first code w the that is now as the 


ys 


e: 28 days” 


Sf, 

‘Strain 


3. 1. —FLEXURAL STRESS “STRAIN CURVES 


somewhat simpler. With allowable stresses 
established, it has for many years yielded safe results, but it i enna a 
in ‘rather unbalanced or non- -uniform 


4: 


"working stress type of of computations. Atl least three factors have been involved —_—- 


First, concrete in tension is quite limited : in capacity. city. In beams, this ten- 


— 
= 
— 
— 
American Concrete Institute. Instead, the apparent simplicity of the straight- 
#4 
— 
Capacity Id larpely 10st l0ading Ib increased to the crac 


ULTIMATE: STRENGTH | 


state. However, with reinforced beams, the cracki cracking load can be as 


_ihighas 50% to 60% of the ultimate load. Indeed, minimum reinforcing is often 
"specified to prevent immediate failure as soon as the first flexural crack ap- 
4 pears. Although ultimate strength design does not eliminate this problem, it 
- focuses attention on ultimate rather than working load conditions; and flexural 
tension in concrete is rarely significant. at ultimate load. 


‘small downtrend beyond the | aa as indicated in Fig. 1. This pattern’ is well 
authenticated, 2 and it is not surprising that any straight line approximation > 


_ ‘Third, shrinkage stresses are often of significant magnitude and at work-— 
ing load are very difficult to assess. Fortunately, when the concrete once 
cracks in flexure the shrinkage stresses largely dissipate insofar as flexural > 
strength computations are concerned. But shrinkage stresses greatly influ-— 
ence the cracking-load itself and, thus, their effect is quite noticeableinre- | 
_ search. Whereas ultimate strength computations d do not | eliminate this | prob-— 
lem, it is less important at ultimate loads. 
Load-Stress Relationships .—Figs. 2 to 5, inclusive, indicate in a general 
- way how widely f, andfg vary with loading. (These are not laboratory data but ia 
rather estimated curves based on several wellestablished check points. x Figs. 
2, 3 and 4 are based on fi. = 3000 psi and intermediate grade steel (fy = = 40,000 
psi), whereas Fig. 5 is more diagramatic. To facilitate comparisons, the axes 
in these figures are given indimensionless units, with the ratio of load to ulti- © 
mate load plotted vertically y and the ratio of unit stress to ultimate stress for i 
concrete (or unit stress to yield stress for steel) plottedhorizontally. There 
— is a major break at first cracking, particularly when percentage of steel - 
small. There is always a second irregularity near — halal wo 
oe In the case of the prestressed concrete beam there is no sharp yield point © 
for the steel, and the steel stress ratio plotted is simply that of f, to the “ 
that exists at beam failure. With prestress, the effects of cracking and of high | 


= _ - eonesete stress tend somewhat to move toget her toward the top of the curve. 


The relationship between loads and stresses is more effectively presented : 


a a ‘working load of roughly half of the. ultimate load develops’ an ‘ that seems to 
‘ be reasonably proportional to load. However, for s small percentage of steel f, 
= far from linearly, and the same is true with prestress 
nce in pointing toward the true ultimate moment. In the case of prestress, 
the dotted extension in Fig. 5 of the working f, curve _to ultimate (in which : 
ae = 1,00) would be overly optimistic. In the case of p = 0.03 (Fig. 4), , the 
dotted extension of the working f. curve would show a very low value; working 
stress methods do yield this low ‘answer, limited by a failure in compression. | 
Thus, to the designer under working stress methods, the low apparent value of 
ultimate load might be increased by adding compression steel. In contrast, the zy 
of actual behavior of this beam would be tofail in tension. Ultimate strength anal-— 


indicates the tension failure and the considerably larger ultimate moment 


InF ig. 7, the variations in fo at working load are wide and give little assist- 
a 


capacity. Hence, the addition of compression steel would add practically no 


om: “Confirmation of Inelastic Stress Distribution in Concrete,” is Eivind Hognestad, a 
“ASCE Paper } No, ST Vol. 83, ‘March, 1957, 
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ULTIMATE — 


shown in Fig. "6. should be noted, however, that whether initial weak- 
a ness is in compression or in tension, the final failure is always a crushing in — 
compression. This is a key point in ultimate strength analysis. The extreme > 
fiber stress is always less at failure than at slightly smaller loads, with the 
fi, stress s existing well inside the beam a and the lower stress corresponding to 
the maximum strain | (Fig. 1) at the extreme fiber. 1 
_ One must conclude that working load and working stress patterns s mask q 
— rather than point out the true ultimate load capacity. There are too many _ 
different shapes the load-stress curve can take. In extreme cases it can even 
point toward added compression steel ina situation in which only tension steel 
s can add strength. This also 30 suggests that the factors of safety under fixed 


Consistent thus al almost demands based on on ultimate strength, 


FACTORS 


— analysis is able to ultimate 
_ rather closely, it introduces a new problem for the designer. Under working : 
~ stress methods, committees from time to time have set up allowable working my, 

stresses that were assumed to give adequate factors of safety. Under ultimate 


strength design the designer must select proper load factors (or factors of 
Some load factors may ultimately be spelled out in codes, as is now the ; 
case in the ACI Building Code . Appendix. 3 However, just a as allowable oe, 
alone have never provided adequate protection in the case of a cantilever 
- "retaining wall (overloads are always considered), in like fashion no rigid set 
- ‘factor ‘must increase » the pan load tobe carried by a ‘beam. but, when the a 
load acts as a stabilizing force on a retaining wall, the load factor should re- : 
7 Much thought is. currently going into the choice of proper load factors, and 
_ there are many variations in ideas presented. Afewof these will be discussed | 
~ briefly because the designer of the future willbe under the necessity of under- 
_ Probability Theory for Load Factors.—A corollary of the need for ‘under-_ 
_ standing load factors is the need for understanding probability theory, at 
in simple form. Consider ahistogram: of compressive cylinder strengths, 4 


that is, a plot of the number o of specimens yielding each strength, as in Fig. 8 
(taken from another publication. 4) If there are a large number of specimens, ‘ 
the ordinates will form a curve that at least 2 approximates a high point at the 
average or mean cylinder | strength and that shows a series of decreasing 
numbers of specimens for - both higher and lower values of strength, This — 
curve is often approximated by the (bell curve 
ro obability theory, as shown in Fig. 8. a = 
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ULTIMATE 
Another convenient concept used in probability theory is the standard devia- : 


in which m is he mean (average) v value | of e. m - fi represents the difference > 


= between a given sample and the mean, and N is the number of samples. Pos- a 
=. , a clearer way to visualize the standard deviation (in the case | of atrue- 


n 
95,45 


6827% 


of Tests 
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10. —DESIGN STRENGTH RECOMMENDED BY 


and the value that leaves below this range approximately 16% of — total 


_ specimens, , as shown shaded in Fig. 8. If one goes 20 f from the mean, the 
i shaded area reduces to . approximately 2.2% of the specimens as 3 shown cross- 
5 hatched in Fig. 8 and10. The shaded area measures the probability of a speci- 
men having a strength less than the value at the limiting ordinate. = = 
Well-controlled concrete will have asmall standard deviation and a 
controlled concrete a large standard deviation. Both conditions are shown in 
‘Fig. 9 ¢ from another publication? , prepared to show the mean cylinder 
yi strengths ‘required to give an f¢ of 3000 psi in n accordance with the ACI Code 
requirement for ultimate strength design. This code requires that not more 
than one cylinder in ten may fall below the specified f¢ . Fig. 9 also shows a | 
further term, called the and noted as V. ‘The: variance is 
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instead of as -arange in psi (for an frequency curve), The larger V becomes, 
the higher must be the average f, to secure a dependable given —, 
a4 Baker's Factors of Safety. —Several years ago A. L. L. Baker, of England, | 
proposed a system of factors of safety that involve both load factor and 
probability | theory. ail Two separate elements: enter into his system, a reliable 
material strength and a separate factor of safety or loadfactor, = | 
‘His discussion on reliable material strength relates primarily to ‘concrete 7 
_ strength. He suggests that the assumed dependable cube strength be taken as 
the mean strength (of a large number of specimens) reduced by two fon one 
deviations, as shown in Fig. 10. (For the normal Gaussian distribution, one 
‘standard deviation deduction would correspond to a 15.9% chance of smaller 

values; two standard deviations to a 2.2% chance of smaller values.) ' There | 
= is then only a 2.2% (or one in forty-five ) chance that lower strengths will 
exist (assuming thet strengths vary according to the symmetrical Gaussian : 

- distribution.) He states that the two-standard -deviation reduction should not _ 4 
exceed 1500 psi “when the most efficient control is used” and 2000 psi for a 
‘medium control. This material is considerably stricter than 


“Results of failure - 1.0 to 


2. Workmanship - 0.5 to 2. 0 
Precast “factory manufactured” - die -0. 5 
3. Load conditions - 1.0 to 2. 
ns High for simple spans and overload possibilities; — 
ss low for load combinations such asl live e load and one wind = 
4 . Importance of member in structure 
_ Beams may use lower value than columns — 
+5, Warning of failure 


4 numbers in the right hand column are suggested by Baker for the frame 
of an apartment building on good ground in an area in which workmanship and - 
inspection are reliable. The column total is divided by 10 and added to unity 


to obtain his factor of safety: 


= 1.0 + DW/10 = = 1, = 2.0 


- ° “The Ultimate-Load Theory Applied to the Design of Reinforced and Prestressed 
Concrete Frames,” by A. L. L. Baker, Concrete Publications, London, 1956. rrealt: 
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ety (or unity a weighted qua — 
ad For his factor of safet shed by adding to unity a v —_— 
Fo to 2.00 These are established by indicated ranges as 
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‘ternational Council for - Building Research, presented a aud that (abstracted - 


elsewhere® had for its opening statement: 


- (initial) cost of construction, the maintenance cost during the useful life, (and) 

: an cost of a hypothetical insurance policy to cover accident risks, odor 
losses, and possible loss of human life , (all) with consideration to. possible 
of public opinion and similar factors of a nature. 


(Terms in parenthesis added by this author. 


‘This report carefully lists the errors in used in design in- 

" cluding: the size of loads, the way they are applied, and the way the struc- 
ture is analyzed by the engineer; the actual strength of the materials used and 
any decrease in strength with time; variations in the quality of workmanship 

_ and the strictness of control; the permissible probability of failure; and the in. 

ductile or brittle failure behavior of the members as designed. 8 
There is a strong position taken in favor of f using a statistical or ww a 
ay approach to safety, (1) reducing | the assumed strength of materials, and 
(2) increasing the assumed working loads. 
For materials, it is interesting to note that this report6 suggests the use of 

7 a  taaie design strength computatedas the mean strength reduced by the factor _ 
—(- - 0) and further lowered by dividing this reduced strength by a factor of 
1.5 to 1.6 for concrete or 1.15 to 1.20 for steel. l. The term o is the standard 4 
strengths so reduced, loads are to be multiplied by a basic factor 
1.25 that, in turn, is multiplied by other factors to represent less than ideal — 

7 conditions. ‘For example, with only average design and construction accuracy 
~ the factor becomes 1.65, with 15% more to be added if there is a very serious 

risk of danger in case of failure, as for atheatre, ,or adam in a highly - ee 

_ ACI Code on Load Factors.—American feeling on on load factors seems to be 
{ 


developing along lines related to the above two plans, at least in the sense that | 
reductions in nominal material strengths | seem necessary (below those used - 
a date). A Also, increases above service loads are obviously essential. Ih addi- 
tion, the original ACI Code Appendix provided different load factors for dead 
"4 and live load. . There is every evidence that the next ACI Code will continue — ; 
this differentiation between dead and live load factors” although probably in 
different fashion from the present code. Some differences should be expected, — ( 
_ in that the present, code uses only one combined factor to cover both material 
_ deficiencies and overload possibilities, whereas the next code will ll attempt to 
-treatthe two separately, 
_ Another important variable to consider is construction accuracy. How should 
the designer allow for a vertical al displacement of moment steel? One- ~quarter 
in difference in a 3-in. ~slab- isa large difference, and in a 24-in. beam is_ 


«Toad Factors Edvardo Torreja, Proceedings No, 55, November, 
p, 567. 


"ability of and (should) achieve a minimum value (total cost) for - 
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ULTIMATE STRENGTH 
— studies of load factors |; and related conditions are now being n mad 


by the ACI Code Committee in an to limit the probability of failure 


ical 


of withmate by Whitney, gave consideration 


to unit strains as a part of his explanation, but not as part of his ‘method. The 
Joint ASCE- ACI Committee introducedthe strain idea into the balanced column 
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FIG, 13. 13. —INTERACTION DIAGRAM, STEEL ALL FOUR SIDES 


In effect, the rectangular stress” block was taken with stress intensity of 0. 85 
4 ft and with deptha = 0.85 ky d (for fi up to 5000 psi). Later experimental work, 
still unpublished, seems to confirm the fact that the depth of stress block can 
be taken as 0.85 ky d, probably | with some reduction for f~ values greater veal 
4000 psi. This gives a most convenient relationship between strains and load. 

‘ Starting with an assumed neutral axis and an ultimate concrete strain of 
+ 0.003 in. per in. (as suggested by the ACI Building Code), it is very easy to 
establish the ultimate load and eccentricity for a column with steel in any ar- 
rangement, as in Fig. 12. One simply evaluates the strains (and thereby the _ 
stresses and loads) on each set of bars and the resultant compression on the a” 
congress. The applied load must be equal and opposite to the resultant of these — 
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FIG. 14.—INTERACTION DIAGRAM FOR TIED COLUMN 
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compared to approximate solutions for the same cases. 

‘Recently the writer hada research test beam that failed in flexural compres- 

4 sion, but that by conventional ultimate strength analysis seemed to be adequate- 

oly designed for more than the attained load. This beam used high strength steel - 

_ (nominal yield strength of 75 ksi) withno sharp yield point and a gradually curv- 

_ ing stress strain curve in the upper stress level. The steel strain was a 

jae of the large stress and larger s still because of the shape of the stress- 

strain cu curve. i Present | design rules are ‘not adapted to this magnitude of yield — 

-- or yield strain. _ (The ACI Code limits fy to not over 60 ksi for beams) — 
The basic strain procedure described ror columns 
beam weakness logically. 


It appears that ultimate strength theory must develop more into a strain 
theory in establishing steel stresses, at the same time taking advantage of the 
simple stress” block concept in} moving from strains to ultimate stresses and 
loadsonconcrete, 
_ In this sense possibly one can talk of load- stress relationships, but this is 
a far cry from the earlier meaning of that phrase. . Ultimate strength design 
ama by its very nature develop at least semi-empirically. _ It bodes well for — 
future use in new situations that it does ; appear to yield to a i logical develop- 
‘meat from strain relationships. ; 
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internal forces, a sifl me modification, but the method is simple and general, 
stants used may need some m 
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STRUCTURAL DIVISION 


aes -DISTRIBUTION CONSTANTS FR FROM | CARDBOARD ANALOGS 

“o> . 


An method of determining moment-distribution and slope- 
deflection constants for beams with variable moments of inertia is proposed, — 


Simple graphics, , especially evolved | for tl the purpose at hand, are used to — 


represent significant parameters of the constants sought, © Template areas are 
evaluated and their centroids are located experimentally. - This is equivalent, Z q 
respectively, to numerical integration and division of one integrated function - 
by another, Thus, the experimental part of the method serves as a mechanical 
- analog, - The m method is simple, r rapid, _ and i sufficiently accurate for office use. use. 7 


7 
BACKGROUND 


us ‘The author has previously presented a semi- experimental ‘method? of de-— 7 


gration, ‘the method is essentially a mechanical analog. 
Note. —Discussion open until June 1, 1961, To extend the closing date one month, a 


written request must be filed with the Executive Secretary, ASCE, This paper is part 
_ of the copyrighted Journal of the Structural Division, Proceedings of the American So- 


ciety of Civil Engineers, Vol. 87, No,ST 1, January, 1961, 0 
Prof, of Civ. Engrg., Manhattan College, New York, 
give 2 “Moment Distribution Constants from Models,” ” by Otakar Ondra Transactions hia’ > 
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_ Perhaps the 2 _— interesting aspect of the me method is that it affords a ‘means: 
of visualizing the effects of moment of inertia variation upon the factors pre- 
“requisite for carrying out moment-distribution or slope-deflection computa- 
tions. . Admittedly, the preparation of the three-dimensional solids and their 
weighing require certain materials, tools, and instruments which are gener- 


=. unavailable in an engineering office. Bog addition to wood, a saw, a weigh- 


a plane. Cardboard and scissors 
are per instead of “wood and s: saw, and weighing is eliminated. Ratios of inte- 
= or algebraic ous are evaluated experimentally by locating comnguets of 
‘aia templates, readily constructed with the aid of graphics, Because 
neither ¢ special instruments nor skills are required, the new methodoffers def- 
inite possibilities as ausefultoolof analysis, 
—The letter symbols adopted for use inthis paper are defined where 
‘they | first appear, in the illustrations or in the text, and are arranged alphabet 
ically, for convenience of reference, in the Appendix. 


moment induced at the fixed endof the member to the moment causing noe 
at the other simply-supported end. In Fig. 1(a), the beam simply supported 
at a and fixed-ended at b is subjected | to a moment t Ma at a, causing moment | 
at b in which Ce is the carry- -over factor fromatob. 
_ Fig. 1(b) shows a prismatic beam whose depth, d, changes from one value 
to another, Fig. 1(c) shows the relative E I- values 1(d) is the 1/E 
‘Area A in Fig. 1(d) is s defined as A= 
it has been shown? that the carry- over er factor a to be ‘computed 


Ja 4 


(statical) moment of area / 


vow 


— 
— 
> & 
form the various required operations, 
purpose of this paper is to propose a much simpler method wherein 

ing a 

a5 
| 

a 
— 

a in which Q, is the first | A taken about the vertical | 2H 

i _ axis through a, and Jg is the second moment of A about the same axis, Fig. 1(d) —- 
designation “second moment” is preferable to “moment of inertia” to avoid — 

Figs, 2(a), 2(b), and 2(c) represent the elevation, top, and pictorial views, 

respectively, of a solid. Its base is the 1/E I-diagram of the stepped beam in 

Fig. 1. The plane upper surface of the solid is inclined at an angle from the 

Quantity Q, used in Eq. 1 may be interpreted either as — 


sT1 MOMENT DISTRIBUTION CONSTANTS | 


—_ is the first ‘moment ~ area A about a in — 2(b), or as 


which is or en the first moment of vaheme V about a in Fig. 2(c). The 


latter may be expressed as V X in which x is the distance of the centroid of 
volume V referred to a. If the specific weight of the solid is uniform, x is also | 
the distance of the resultant ; body- force, W, that is, the distance of the center 


— 
Using the meaning of Qa and Jaa as s defined 7 Eqs. 3 and 5, the ratio Qa Ny 
in Eq. 1 is equivalent to 1/x. By fashioning the solid shown in Fig. 2(c)- a 


its appropriate weighing as has been described, 2 - distance x may be eaeneinnen 
and > Eq. 1. 


3(a) shows factors I and k: x of the in Eq. 3 juxtaposed 
inone plane. If the } product of each pa pair of colinear ordinates is effected, (1/E I) 7 
(kx), laid off in Fig. 3(b) as ordinate and connected by a line, the plane area — 
bounded by the latter and the axis a-b is equivalent to the volume V in Fig. 2 (c) 
and in Eqs. 3 and 5. it follows that the body- force distance, x, in Fig. 2(c) ll | 
the same as the centroidal distance of the area in Fig. 3(b), _ 7 
‘The distance x may be determined experimentally by cardboa. d 


complished by balancing the template on a knife-edge, or more accurately by 
‘suspending it on two non-colinear strings as in the familiar experiment of el-— 
ementary mechanics. The transformation of the solid into a two- dimensional 


a 
3) 
— 
"which is (1/k)-times the volume, V, of 
— 
a 
ia 
— 
— 


APS 
the relation between the carry-over factor, Cap, and the 
metric parameters L and X of area V in Fig. 3(b) may be derived more readily 


ed direct methods. For example, in an earlier paper of the author’s, 2 previ- 
ously referred to, his Fig. 1 external work 


b) 


‘Se 


paper is the load on the eiiltesi beam for the real — in the . earlier Fig. _ 

1(d). It follows that the ratio of the conjugate- beam reactions, and hence the — 

location of the. the centroid of V, determine the carry-over factor, 
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-MOMENT- DIST RIBUTION ow constants 


‘order to m make the cardboard template, both tasks may be slbeeaiad ot 
rapidly accomplished in a single graphic operation, 
Graphic multiplication is based on the geometry of similar triangl in 
4, x is the axis of and functions f(x) y=kx=x are 


shown, ‘To m multiply £(x2) (x2), ‘set 


in which the divisor on the left, 1, is a unit of x. If (x9) is is projected horizon- 
— tally into 2 on the vertical line x =1, the ray drawn through the origin, 0, and 
and extended to the vertical at x2 determines (1/k)P9=2p 9. 
In Fig. 5 the graphic multiplication is carried out for eight subsequent. ab- 
scissae x. The manner in which the line P = f(x) x = F(x) is obtained is - 
- lustrated for x-values of 2 and 5. The ordinates of P are drawn full- scale, that 
_ is, they are read on the same scale as was used to plot the ordinates f(x). 
ae the scale of the P-ordinates may be controlled as desired, for bya] 


ample, if P is to be obtained to the one- quarter scale, Eq. 7 may be written in 


> 
a = 4 and rays connecting the the extended, are drawn 
shown for 3 8. The ordinates of the line so constructed, 
ra == 1 ), are ‘the p- ordinates to the one- -quarter r scale, 


To multiply f f(x) (kx),a ‘comparison of Eqs. 7, 8, and 9 shows that k may 


be interpreted as a scale factor of the P- ordinates, ‘However, it is significant _ 
to note that the abscissa of the centroid, x, of the area under the P-line is in- — 
dependent of the value of k ar k and any other scale factor ( Eq. 6). 


herein considered the scale factor need not be taken into account. “In practice 

7 : the choice of the scale is governed by the desirable graphic magnitude of the © 
maximum P-ordinate. In general, will occur at some fairly high x-value 
coupled with a large value of 1/EI. Acursory graphic multiplication for one 

4 or two suchtentative x-values will ascertain the scale fa factor necessary i in or- . 
der to confine the plot to a desirable size. me 
The graphic procedure may be further simplified by noting: that it is actu- 
ally unnecessary to plot the line 1/EI=f(x). Values of f(x) may be plotted 
pen on the vertical ae 1 (or r any other vertical if a scale factor is —— 


The variation of 1/E 1/1 values, absolute, or relative, fundamentally 
eachof the moment-distribution constants. In standard methods of anal-_ 
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before, 
the latter and the P-values are found as 
— 
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‘January, 


(4. 


= 
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sis, values of I and their vecinvecala, 1/1, must be computed, then tabulated 


r plotted, and several properties of the analytic function 1 fas f(x) or its geo- 
metric diagram are evaluated. , It will be shown subsequently that for the pur- 


graphic of the ‘cardhoard template. 
the possibility of numerical 
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MOMENT-D -DISTRIBUTION CONSTANTS a 


= 4(x), 8 


rewrite in the form ‘a 


“4 In Fig. 6, bisect the aaa angle between coordinate axes x and I, and draw 

= the horizontal line I= 1. The construction of Eq. 11 is illustrated for values © 
Xg = 2 and Ib = $(xg). Project (x9) horizontally into 2 located on the bisec- 


tor, hence vertically into [2] located on the line I = 1 . Through the origin, 0, _ 
of the coordinate axes, draw ray 0- . and extend the same to its intersection — 


which is evident i in satisfies 


the I- ordinates into abscissae of the horizontal axis I=1, 
quently, the I-diagram and the need not be drawn, yoy I-values cor-_ 


The rest of the construction remains the same as in Fig. + analy — 
simplified graphic procedure is illustrated in Fig. in which a beam 


mt with a parabolic intrados and a constant - depth haunch is chosen 1 a8 an ex mangle. 


: 
relative I- values were computed for sections a unit distance 
and recorded in the figure. A suitable scale is chosen forthe I-values and the 
horizontal I-axis is drawn at the distance of I = 1 below origin 0, Values of I 
are laid off on that axis and labeled [0] , 1], [2], and [3], ete., thus indica- _ 
ting x-values of 0, 1, 2, and 3, etc., with which they are associated, Rays drawn - 
the origin, 0, through points [1], [2], and [3], etc., intersect verticals at 
abscissae x = 1, 2, and 3, etc., at points r which determine the ratios x/I= = a Ty, AY 
2/19, and 3/13, etc. The segmented line connecting the thirteen points obtained 7 
int the foregoing | manner is the line ‘R = = dotted | of 


_ yon the one-half scale, either the ordinate I « = 1 is halved andthe I- -scale 
_ kept as before, or the I-scale is doubled and the ordinate of the I-axis, T= i 


ta 

% 4 = _ To divide values of x by — 06a 
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MOMENT -DIST RIBUTION CONSTANTS 
graphic con: 
should be done on a sheet of condensed: bristol board, or similar paper of sub- — 
‘stantial but uniform thickness. The required template - is obtained by cutting» 
out the area bounded by the R-line and the x-axis, The centroidof the area can © 
accurately and located inthe manner ‘shown schematically in Fig. 
: ‘The template is hung on a pin or thumb tack, 1, passing freely through ae 
es and embedded in a smooth vertical surface. A weighted thread is oe ae 


ped around the p pin, Care must be taken to avoid friction between the template 
and the pin, the vertical surface behind it, and the suspended thread, This ca 
a checked by slightly displacing the suspended template to one side of its po- 
sition of em, and when released, the template should oscillate freely. 


plate is marked by a dot uinice along the thread, and the line 1-1 is ——- va 

‘connecting the dot and the pin hole. The same procedure is repeated using pin ‘e 

_ 2 placed at some different position with respect to the template and line 2-2 is 

_ obtained, _ Evidently, the point of intersection of lines 1-1 and 2-2 locates the — 
center of gravity of the template. The accuracy may be with the 


bs a third pin, 3, as is gga in the third diagram of Fig. 8 . 


"gravity ‘and thickness, the center. of gravity also locates the ¢ centroid of | 

area. The centroidal distance, x, is scaled off | and its value substituted into | 1 

Eq. 6. In the example under consideration, by actual measurement, x was found — ; 


to be 0.555, so that Eq. 6 


By equivalent analytic computation based on subdividing the 1/1- -diagram 


into twelve strips of equal lengths Ax and computing Zz Ax ‘I and « Ja 


> 
0.8 81 


ete., the templ: required 45 min of time. ‘The 
- merical counterpart was accomplished in slightly less than 1 hr, The task of © 
7 _ Computing the I-values, , which is common to both methods, was not timed, The 
: analytic solution is approximate to the extent that finite differences rather than 
_ . integration were used. The former approach is the standard method of deter- — 
_ mining carry-over factors for beams with continuously varying I. Since twelve _ 
~ strips were also used in the analog method, the two- solutions have a eieenal 
basis for the comparison of accuracy and time expended. Pee 


In general, aside from the small errors inherent to any graphic or experi- 


mental method of analysis, the accuracy of the analog method depends some-_ 


what on the number of segments Ax used in eng the » or R-lines, As 
this number is the area under either line, (x Ax) ‘approaches 


(x 4 dx) It will, be actual use of the method 


x 
| 
a 
as 
q 
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i _It is worthwhile to note that the difference between the analytic and analog- [iim 
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MOMENT- DISTRIBUTION CONSTANTS 
that the graphic nof the template i is so and rapid that it 
= whether 10, 12, and 16 or more elements are used, Since the work of . : 
balancing the template is the same regardless of of its make up, the degree of ‘S 
the controllable part of the accuracy depends on one’s ;decision as tohow — 
elements will be used, | Obviously, this number should be larger f for beams 


eT 


di agram 


44 


in = @ 


witha a rapidly changing depth. In practice, the task of computing I- values ues at 
closely spaced sections may be simplified by constructing the P- or R-line for 
larger intervals. ‘The line is then —' by a continuous curve drawn tan- 7 
to the segmented ling 


— 
— 
. 
im 
had 
— 


y, 1961 
anu ry, 


"small e errors in the graphic experimental nature of the analog 
The factor, | S, at the simply supported end ofa | member 
 % defined as the moment required to rotate that end through a unit angle while ‘ 
_ the other end is held fixed. In Fig. 9(a), Mg is is Sap, th the stiffness factor at end 
a of the beam a-b, when @, is one radian. 
_ Inthe Hardy Cross method of moment distribution, 3 an unbalanced moment ha 
acting ona joint is distributed among the members meeting at the joint in pro-— 
portion to the distribution factors, S,/ZS, in . which Sn is the stiffness factor at 
the joint-endof the member under consideration, and ZS isthe sum of the stiff- 
ness factors at the joint-end of the members framed into the joint. | 
~ It has been shown2 that the stiffness factor — be ee from the for- 
mula 


(14) 


“onl 
“The te term in ‘Eq. 14 may be as Qa (Qa/Jag), or, 
_ Eq. 6 , Q,/%. But Qa is the first moment of area A about a (Eq. 2), and may 
therefore be set equal to A Xc, in which x¢ is the distance of the centroid of 
area A from a. Consequently, Eq. 14 may be written in in the form 


oF Distance x X appeared in Eq. 6 and its analog- evaluation has already been ex- 


- plained and illustrated. Quantities A ar A and dX¢ may be be found in the e following m; man- 
‘Fig. 9(b) shows a W beam withcover- -plate CP, tie-plate TP, and bearing-_ 
plate BP. Relative I- and 1/I- [-diagrams are shown in Figs. 9(c) and 9(d), re- 
_ Refer to the schematic diagrams in Fig. 10, Beginning with a nates 
sheet of cardboard, a~b-c-d, of area B square inches (say, 10,00 in by 8.00 in), | 
Fig. 10(a), the 1/I- -diagram is drawn to scale as is shown. Next, the template 
for area A is cut out with a razor blade. _ The centroid of A is located in the 
same manner as was one in Fig. 8, and the centroidal distance, Xe, Fig. 9(d), 
A Likewise, the — of template e-b-c-d of area (B - A) is found experi- 
_ mentally and the centroidal axis is drawn as is shown in Fig. 10(b) (dashed 
line f-g). Next, the template of Fig. 10(b) is cut into two parts, Figs. 10(c) _ 


_ and 10(d), and the areas of the parts are denoted by AL and Ag. The centroidal - 


distance of Aj, x4, is found by direct measurement; that of Ag, denoted by ry XQ, 


na: “Analysis of Continuous Frames by Distributing Fixed- ‘End Moments,” 
Cross, Transactions, 96,1932,p.1. 


- 
he 
— 
— 
— 
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from which 


- Areas B and Ay are rectangular and can be computed from measured sides 
with a high degree of accuracy. Centroidal distance i) is obtained with pre- 


y direct measurement. Centroidal distance xg can be e determined very 


cision b 


accurately b by the method of Fig. 8. 
Aside from purely numerical computations, area A could also be determined 


_inother ways (for example, from the measured size and weight o of sheet a- b-c-d = 
and the weight of template A). However, aside from being impractical, ‘it was 
conclusively established that the weighing is frequently unreliable. 
_ i. it is advantageous to determine the absolute rather than relative value of MN 
“area A if fixed-end moments are to be found later. Fixed-end moments re- 


sult from loads placed onthe member consideredor from joint- displacements 7 
caused by settlement or temperature change, 


Relative values of areas A for two consecutive members free of loads and 
_a transverse joint-displacements, such as members c- b and c-d in Fig. 11(a ), 
suffice for the moment-distribution method. 
‘To determine relative stiffness factors of members c-b and c-d, lay out = 
‘their I-templates to scale as is shown in Fig, 11(b). (A simple case of 
I-variation is chosen to illustrate the procedure, Obviously, an equivalent nu- 


merical solution would be more here. ) the c centroid of (Acb + 


Separate areas Ach and i i cutting through c, locate their centroids by the 
method of Fig. 8, and scale off centroidal distances Xeb and Xed- Taking first 


moments about c yields 


_ CARRY-OVER AND STIFFNESS FACTORS AT END b 


a MOMENT -DISTRIBUTION CONSTANTS BT 
— (A A) ee (16) 
— 
&g 
— 
— 
aa 
— 
&g 
Applying Eq. 17, relative stiffnes 
— The fac .6and14, 


4 


| 

— 


the member (its I- values) i is symmetrical shout the 

_ mid-span, the carry-over and stiffness factors are the same, reapectively, a *- a 

ends a and b, If the member is non-symmetrical, factors Cpa and Sba may be © 
evaluated by applying | to end bthe procedure outlined for end a, that is, the 
elite of abscissae x is placed at b, the corresponding P- or R-lines are plot- — 
ted, etc, It is emphasized that t the centroidal distance, x- x= barred, when n refer- 
red to b, is ot equal to L-£ x. 
factors and may be determined 1 more from the 


- Substituting for Je from Fa. 26 Eq. 27 


Combining Eqs. 6 and leads to 


“Substituting i into Eq. 22 fo for py and Jp given by Eqs. 25 and 30 yields 


= 


and 
— 
In Fig. 12, c-c centroidal axis of area A, parallel to a-a and b-b. 
for J, expressed by Eq. 29 into Eq. 28 results in 


(X ; : 
; 


MODIFIED STIFFNESS FACTOR 


The modified _——— factor at end aof a simply supported | beam a-b, Stab 
is the moment which, when applied at a will cause that end to rotate through a 
unit angle, The use of the modified stiffness factor expedites moment-distri- 


computations for structures waich are simply supported. 


"computed from the formula 


| 


Sa 


settente, the formula for the modified stiffness f factor at end b is 


en — 


‘SHEAR STIFFNESS AND FLEXIBILITY 
The shear stiffness of a member a-b, Sg, is 5 the transverse force | Pp required 
to causea unit transverse displacement, A= 1, between ends a and b, heldfixed 
against rotation. In Fig. 13, force P is the shear stiffness when A is unity. _ 
It has been shown? that the by the formula 


aT 


Substituting for ‘obtained 


» Flexibility, F, is tt the of eer ‘stiffness, that is s, the transv erse 
Geplacement, 4, , between ends a and b caused by: a unit transverse force 


Spa = — 
— 
le 
— 
4 &g 
e — 
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7 
Relative shear- stiffness values may be used in sidesway anal- 
of rigid frames, . Absolute values must be ‘used when 
for the effect of settlement of supports or change. 


_FIXED- END IN 
Concentrated Load.— Column-analogy formulas4 will be 
cardboard-analog for -~end moments, In Fig. 14(a), 


‘eees area under the (M/E I)- diagram for beam a- -b cut back to a Pe 
_ determinate one by removing the fixed support at aj it is also the volume oo 


the solid, Fig. 14(b), whose base has widths equal to the (1/E I)- -values and oe 


> — upper surface rising as Px; _ 
= KL + x, 


is the y of W, that that is, of the centroid 


eccentric 


Pons values of 1/E I * values of x. The graphic multiplication is performed - 


in exactly the same manner as was done in Fig. 5 except that it begins at point 
_k. In Fig. 14(c), two W-templates are constructed for load P = 1 — at ky 
- (solid line) and at kg (dashed line). If preferred, templates W, and Wo may 

also be obtained by graphic division of x-values by) the corresponding a 
of EI. In either case, relative E I- or I-values 1 may beused, “_ 

_ The W-template is cut out, its centroid is located by the method of a 8, 
is scaled off. Assuming that is known from 


- centroidal distance Xw 
stiffness-factor eccentricity, e, may be readily computed from 
Be . Grinter, The Macmillan Co., New 


on 
4*Theory of Modern Steel Structures,” by L The 
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7 The function W may we ‘evaluated by the ‘method of Fig, 10 a and Eq. : 
which W replaces A, 
Substituting into” Eq. 40 for e given by Eq. 43 and for Je obtained by com 
bining Eqs, 26 and 29 leads 


Funetion x/E I) in Eq. 42 may be to any suitable 
scale that will yield a W-template of desired proportions. It will be recalled 
a that in graphic multiplication, the scale of the product - function is governed by 
the selection of the abscissa, x=1, Fig. 5 or Fig. 14,x=4, Fig. 5, etc., onto | 
_ whose ordinate the dependent variable, 1/E I, is projected. Similarly, ‘when 
graphic division is used, the choice of the horizontal line, I = 1, Fig. 7, I = 1/4, 
UEC. along which the I-values are laid off ee the scale to which ordi- Pi 
_- The value of W as used in Eq. 44 is the — obtained experimentally by the 
_methodof Fig. 10, multiplied by the scale factor (for example, 4, if the (x/E I)- 
- function is constructed to the one- quarter scale), and by the load intensity, J 
>. _ The term W, defined by Eq. 42, is a function of 1/E I and L, as is A in Eq. 
4 44. It follows that the scales for 1/E I and (L - KL) in the W-template [ Fig. . 
— *14(e c)] should be the e same, respectively, as the scales for 1/E I and L in the 
: _A-template, Fig. 10(a). If different scales are chosen in constructing the two | 
—— the scales must be correlated when evaluating Eq. 44. _ The corre-_ 
lation is effected by applying an appropriate scale factor in Eq. 
Suppose that in the case of template A, the span, L, of the beam is laid off 
as \a-inches, a unit of 1/E I (absolute or relative) is plotted as iq-inches, 
4g and that in the case of the W-template, &- K L) is laid off as w- inches, and — 
a unit of 1/E I is plotted as iw-inches, | Then the numerical values of A and W 7 
y as found independently of each other in the manner described in Fig. 10 ge f ; 
converted to the true ratio- value, | W/Al, which be used in Eq. 44 as 


re 
ow — 
wAw 


An application of Eq. 45 is illustrated in the following 
“gake of simplicity, assume I= 1/5 = constant, and K = 0, 5, Fig. 15(a). Let * 

= 10, i, = 0.4, Fig. 15(b); Aw = 8, and iw = 0.3, Fig. 15(c), 
: ‘The experimental evaluation of A by means of Eq. 17 and the procedure il- — 

lustrated in Fig. 10, and of W determined in the like manner, would yield A and 


Ww expressed in square inches. In the simple case under consideration, how- — 


— 
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_ 4)(10 A, - 0.52 4g 
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laid off is drawn at I=1, ‘This becomes | evident when one realizes that, ae 

doubling the I-scale would move the horizontal line to twice its present distance 

it below the x-axis in Fig. 7, as it would also move t the bracketed points plotted = ‘ 
= that line to twice their present distance from the R-axis. Consequently, | 
_ if graphic division is used to construct the W-template, Fig. 15(d), the hori-— 
zontal line I = 1 ane 1 and tl the segment I= 1/5 t [0 to may be laid id off 


| 
a or _ Reference to Fig. 7 sho edbythe scale - 7 
chosen for the I-values as long as the horizontal lineon which the I-values are 
— 
— 


MOMENT - DISTRIBUTION ‘CONSTANTS 
to any scale, template is affected only by the Aw-scale, 
iw in Eq. 45 and the subsequent computations is replaced by 1. Based on the 
alternate W-template in Fig. 15(d) and simple computations, W is 160, and 
| is found to be (1/8) ) PAw as before, 
experimental practice, however, the W-template in Fig. 15(d) would 
be of impractical proportions ; and therefore, a scale reduction would be desir- 
my able, Using a scale-reductionfactor of 3/10 (dashed construction lines in =: 
15(d)) results in a template which is identical with that of Fig. 15(c). al 


value to be substituted sis WwW on the right side of Eq. 6 would be 


he value of iw int the same equation 


Refer againto Eq. 44, _ Although the second fractional expression on : the right 7 
side is dimensionless, consideration must be given to the x-scales (A-lengths) 2 
; of the three templates, A, V and W from which ‘the four parameters were deter- | 
mined. In general, any three parameters in the fractional expression under i 
consideration may be adjusted to the scale of the fourth, or all four parameters a4 
may = converted to values corresponding to unit span-lengths, A =1, The _ 
span- -lengths are the lengths of tl the respective templates with the exception of 
for which the is the template length increased by 


-scales have no effect upon the position o: of the centroids in the 


_-X=direction so that their correlation is not necessary. 


_ Uniformly Distributed Load, w, over Entire Span. — The fixed-end ‘moments 


caused by a uniforml distributed load, Fig. 16, are 


Ww 
and 


=~ 


‘st atically condition nah removing the ‘Support | at a 


Eq. 48 may be writtenas 


— 
> 
a 
— 
> 
q _ Substituting for J, from Eq. 29 — 


x 


= 
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- Substituting for Ww from Eq. 50 in Eq. 46 and replacing Ie by its cuivaent ob- + 
‘tained by combining Eqs. 26 and 29 leads: 


Sesematees X and x, have already been encountered in the analog-evaluation | 
7 of the carry-over and stiffness factors, and will therefore be considered known, 
_ The only parameter that remains to be determined is e, the eccentricity, which — P 
is defined as the x-distance of the | centroid of the (M/E I) I). ba gy _— the 
According to Eq. 48, ‘if a template representing ‘the integral 
prepared, the x-distance of its centroid from any vertical reference axis is the 
same as the x-distance of the centroid of the (M/E I)-diagram | from that axis, 
‘The template can be constructed from the V- -template (Figs. 3, s and 7) — 
multiplying its ordinates by the corresponding abscissae. This multiplication 
by the x-values amounts to considering the V V-template area to be what ar area A 
_ The graphic construction for the beam and V- -template of Fig. 7 is shown in 
Fig. 17, in which the ordinates are obtained to ‘scale. Accordingly, 


a 


"locating its centroid the ‘experimental method of Fig. 8. centroidal 
_ tance, Xw, is scaled off. Since this distance is independent of factors 2, E, 6, 
a and w, each of which pertains to the ordinates of either Wor W (Eq. 52), dis- 


tance XW is the same for both WwW and W. Therefore, _ eccentricity, e, of w 


which, ‘substituted into ™. 51 results in 


X- 
= 


"Although it it i is advantageous to utilize the V- template. in the foregoing con- 


struction, atwo-foldgraphic multiplication by x may be readily performed with- 
out actually drawing or tracing the V- outline. ‘The procedure is illustrated - 

_ Fig. 18 in which f(xg) is multiplied by x = 2 once, and the product, 2-p'9, is 
thereupon multiplied by x = 2 once more, , 80 that 22 {(x9) = 2 - p's. "Points 
p" ‘sO obtained and connected define the boundary of the required template, AD 
scale 1 reduction may be effected in either one or both multiplications by pro- 
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Couple or Longitudinal Loa ~The fi fixed- moments caused by 


ii couple or moment, C , Fig. 19, are 
Ww 


55) 


in which 


the moment in the cut back toa statically determinate 
one is constant, the eccentricity, e q e, of 7% the distance of the centroidof area 
_A' from axis c-c, Centroidal distance x’ b for template A' is found experi- _ 
mentally by applying t the peocedere of ‘Fig. 8. 8. ‘Then 


_ Relative Transverse Displacement between Ends a and b, Fig. pore —It can y 
be readily shown (for example, by conjugate- beam or column- -analogy methods 
= analysis) that fixed- end r moments caused by transverse displacement, A, 


ACL 
My - 


Cx 
= 
— 
| 
— 
— 
Eqs. 57 and 58 into Eq. 55 leads to 
— 
— 


i 


The moment required to hold wiht b fixed is 


In general, the values of t the constants determined: with the ; aid of cardboard — 2 
templates must be adjusted to express the constants of the prototype. The fun-— 
damental dimensions in the construction of templates are lengths, Lor 7 
VE I-values, plotted as scaled lengths in a rectangular coordinate system, 
Centroidal x -distances, scaled off from the be 


Stiffness Factor r 


Modified Stiffn 


Shear Stiffness 


g? 


3 


re 


Fixed- End wload 

Fixed-End Moment, C-load 

Fixed Fixed- End Moment, A-displacement ) 
Fixed: d-End Moment, 6, ,~displacement 


in the prototype. Since ratios of like quantities are re dimensionless, they retain 
the same value when applied to the prototype. Areas of templates, computed ; 
_ in square inches, must be multiplied by the scales of the two fundamental di- 
Let 1 in, of template length r represent m- inches of sp: span n of the prototype, 
_ and 1 in. of ordinate represent 1- (1/pound- ound-square inch inches) of VEL ‘Table 1 i _ 
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"MOMENT DISTRIBUTION ( CONSTANTS 


the various moment-distribution constants determined by “cardboard 


“analogs, and constants adjusted for the prototype. 


-CONCL LUSIONS 


wa bya group of fifty senior civil engineering eames at Manhattan 1College, 
New York, N. Y. The moment-distribution constants for one of two ee 
‘structures (a three-span continuous steel beam with fifteen different I-values, | a 
and a two-span concrete bridge frame with varying I) were determined by each : 
Student with the aid of cardboard templates. To provide a check, several stu- 
dents and the writer also ‘computed the constants by analytic methods, _ ah 
_ In general, an excellent correlation was found to exist between the experi- 

mental and analytic results, Occasional when were 
traced to one or more of the following — - 


"Gross errors in the initial computations of I- values, 
of unsuitable template material (corrugated ordinary 


3. _ Choice of very small m-scales (templates less than 4 in, oa” 
4 Poor drawing technique, crudely | constructed or cut t out templates. em 


the aforementioned sources the correlation be- 
tween template and analytic results was generally well within 5%, with the ex- 
ception of fixed- end moments caused | by c concentrated or distributed loads where | 
differences of uy up to 10% were not uncommon, In the latter case, the discrep- 
ancy is not necessarily or entirely due to the analog technique, An inspection = 
of the column- 1-analogy formulas on which the working formulas of the card- 

- board analog ‘method are based, shows that differences of quantities. of com- 
parable order of magnitude are involved, 2 Thus, both the analytic and experi- 
mental results are sensitive, respectively, | to the accuracy ° with which the sev- 
_ For all practical purposes, the cardboard analog- method is sufficiently ac- 

curate, It is also highly ‘reliable s since little opportunity e exists for gross 


‘beams with constant or variable moment of inertia, and of statically deter mi- 


or or indeterminate st structures in sin general, 
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=t “The Determination n of Moment Distribution Constants. of Members w with a 
Variable Moment of Inertia,” by °. Ondra, thesis presented to Lehigh Univer- 
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PPENDIX, — NOTATION 
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SOLUTION OF BEAMS WITH VARIABLE MOMENT OF 


L, PEI,* M. ASCE. - —Many interesting questions were up by 
discussers. _ The following is an attempt to clarify these points: 
_ An Alternate Derivation of the Matrix Formula.—Instead of the finite-dif- 
ee approach, the problem may be considered from the integral-equation ‘ 
approach, The ‘equations of the simple beam are two second~ — 


ax’ 


Since | een two equations are or in form, they y can be solved inthe | same 


manner in succession, as two integral equations 


= fale, z) q(z) dz .. 
| 
where A A(x, 2) is known as the ais Function of the differential equation 
he 6). The problem is. then reduced to the search of the Green’s Function, 
The matrix [A] is, in fact, this Green’s Function in finite-difference form, 
and actually represents the exact solution, in closed form, of the finite-differ-_ 


sert a weighting function, or weighting matrix [ w] (called concentration ma- 
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Szabo? Ss weighting matrix depends on the values of ‘moments of inertia, and can- 


not be combined with [ A] before the writing of the machine program. Further-_ aa 
more, the values of moments of inertia at half-intervals are required, which | 

is awkward, Finally, this weighting matrix displays a peculiar asymmetry, © 
which should not be present in weighting matrices, 
Weighting matrices (Eqs. 32 37) should give better answers than 


~The increased accuracy is obtained at the = 


is to be done on the computer. . To be consistent, “it would seem that the same 
weighting matrix should be used both in Eq. 30 and Eq. 31. Dis 
ea Computer Program.—Fenves raises the question of how many storage lo- 
am as required to progr program the writer’s formulas, It appears at first that | > 7 
en n2 locations are needed to store [ A]. Actually, it is not necessary to 
[ A] which can be generated, as it is needed, by the 


< 


if) 


= developed by the University of Michigan Computing pee for the 7 
= 704 computer which is shown in ~ 6. Such a program can be ew 


It oad require considerable sophistication to write as short a ‘program using — 

(New mark's Method. —Newmark’s Methodgives an 1 extremely efficient algo- 


— 
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ere locations, The program is quite short, and did not use special matrix subrou- — 4 a 
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trix, algorithm is described recently by (17). The writer’ s scheme 
“plugs given values into the formulas 


In spite of this, Newmark Ss Method cz can 1 be | faster if all the deflections are re- 


BEAM DEFLECTIONS BY MATRIX FORMULA 


READ FORMAT INPUT 
READ FORMAT INPUT1sQ(1)eeeQ(N-1) 
PRINT FORMAT TITLE 3 
NE 


EQUIVALENCE 


“THROUGH ABLE» FOR J 
END_OF CONDITTONAL 
THROUGH BAKER*FOR 


BAKER» S=1919Se Ge (N= 1) 


WHENEVER R 

YER) ZY(R)+R# (N= 5*ME) 


— a END OF CONDITIONAL 
DIMENSION) ¥(100) 


“PRINT FORMAT DATAsY(I}eeeY(N=1) 

VECTOR VALUES TITLE=$1H1 S38 »44HOEFLECTIONS OF SIMPLE BEAM BY 
“VECTOR VALUES INPUT =$S409F10. 
VALUES TNPUT1 = 

END F PR RAM 

OF PROG 
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Clough" Met 
garding his method (16). . His derivation of flexibility matrix is certainly not 
as straight forward as he e claims. . After obtaining a 2n by 2n flexibility ma- 

- trix, he then arbitrarily reduces it to an n - n matrix for computational use, 


hod.— The writer cannot: agree with -Clough’s comments 
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DISCUSSION 


Remarks .—The writer has not been able to locate Szabo’s 


(14), His later papers (12,13) deal with the analysis of orthogonal — 
"structures, and do not seem to discuss deflections of beams | with variable mo-— 
a Reggini kindly called to the writer’s attention two papers (10,11) published _ 
in Argentina where the same idea is explored. The writer wishes to thank ot 
the discussors for their interest in the paper, 

Acknowledgments. .—The computer program mentioned herein was written 
while the writer was serving as Visiting Associate Professor at the University | 
of Michigan under the sponsorship of the Ford Foundation Project on Compu- 
ters, Donald a Katz, Director. - Carnahan gave assistance in MAD program- 7 
ming techniques, Their generous supports are gratefully acknowledged, 
_ Errata.—Page 11, line 1, Eq. 23 should read Eq. 22, Page 11, 5, 
the reciprocal of [h/E is given instead of [h/EI]. 
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poy to offer a general reply. Many ‘tests have shown that full ultimate re- 
"sistance can develop at all critical sections of a statically indeterminate ro 
_ inforced concrete structure before failure occurs at one of them, But compu- 

: tations and tests also show that cases are possiblein which some critical sec- 
tions fail by excessive strain, before others have attained full strength. Hence, 
- critical sections must be designed to be able to deform inelastically before 
‘failure, an amount compatible with the relative stiffnesses and deformations 
_ of the influential frame members. However, the (n + 1)th hinge by appropriate 
_ reinforcement design can be made to form at a load only very slightly tages 
_ than that which produces the nth hinge. Amore difficult problem for the de- 
signer isto decide what alternative distributions of load may occur and be crit- : 

Tieng Here the statistical approach is of great interest, and could be used as 
_ an aid in deciding what t loading conditions and hinge positions to adopt as sound 
design criteria. Having decided on certain load distributions as possible crit- 
ical, each case requires investigation to ensure adequate strength whenthe nth 
hinge | forms, It would not be wise to risk the possibility of failure due t i 

_ incompatibility. Research m may, of course, , eventually show that the 
_ recommendations of design criteria are a little too cautious, and that further 


vember, 1959, by A. L. L. Baker, 
, Concrete Tech., Univ. of London, London, England, 7 


osure by A. L, L. Baker 

L, L, BAKER,+—The contributions of the several discussers are greatly 

appreciated. Since the discussions refer to design criteria it may be appro- 
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Closure by Ardis White 
ARDIS WHITE, 1M. ASCE.—The “discussion by King has considerably 
creased the value of the original paper, as he added a number of ideas which — 
were very pertinent to the subject matter, but were not included in the paper, | 
King has asked about the time and equipment involved in preparing data for 
7 entry into the computer, A special form is used for presenting the data, and 
: about 5 min is required by the engineer to fill out this form, An additional 10 
min 0 or so is used by the key-punch operator in punching data cards for the 
IBM 650, Feeding the program deck and datacards into the computer requires 
‘ about 5 minor less, and the actual computation consumes about 10 ‘min, '_Print- 
ing out the results requires about 5 minor less. __ 
The program is coded in SOAP, including SIR, with sine, cosine, and square | 
— root fur functions, About 12 man-weeks were involved in writing the program; 
coding : and debugging each required about 1 additional week, epee: 


4 pe. The writer wishes to express his appreciation to Arthur McMorris, eoodl 


4 November, 1959, by Ardis White aaa 
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NTS WITH NONPRISMATIC MEMBERS?» 


L, LEE,!M, ASCE and F. P, WIESINGER, 
with Rogers — that approximate solutions are often necessary and sometimes 
q indespensable i in engineering practice. The approximate method used Py — 


as an illustrationof his remarks was originally suggested by Saliger.3 In reine 
~ method, however, the location of the imaginary hinges in the chord mem- a 
J ber (or cross member) should not be in the ratio of the squares of the ett 
= esses of the adjacent cross members | (or chord |members), as § stated in Rogers’ 
discussion, but in the ratio of their square roots, The latter was a 
z used in the numerical workof the discussion with the exception of the hinge die 
cation in the cross members, which cannot be verified. Sareea eo eG 
“a It is most unfortunate that the method | presented inthe paper under r discus- 
.- which deals only with Vierendeel bents with equal chord stiffness and 
symmetrical cross members, is not applicable to Rogers’ example i in which 
_ the chord stiffnesses are different. The latter, however, can be readily treated 
by means of another method of analysis suggested elsewhere by the writers. 4 
_ It should be emphasized here that the method presented in this paper involves — 


‘simple expressions and yields correct solutions with a speed of | convergence 


rapid enough for office practice, 
Corrections. —In the e original paper, on the ‘second d equation from the 


tions of all the forces ‘should be reversed, “ id 


December, 1959, by L, Lee and F, P, Wiesinger. 


prof, of Civ. Engrg. ‘Northwestern Univ., Evanston, mn. 7 
Natl. Science Foundation Faculty Fellow, ‘Northwestern Univ,, , and Asst. Prof. 
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“Der Eisenbetonbau,” by R, Saliger, 6 ‘Aufl., Leipzig, 1933, also Beton- 
-- Kalendar, 1945-1950, XXXIX Jahrgang, Wilhelm Ernst & Sohn, Berlin, 1949, p. 237, 
> -* “Moment - Distribution Analysis of Vierendeel Bents with Inclined Chords of Un- 
equal Stiffnesses,” by S. L. Lee and F, P, Wiesinger, Proceedings, I.C.E., Paper No, 
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-R, S, ALWAR,*—A interesting iterative procedure for the evaluation of 
stresses in beams on elastic foundations with discontinuous loadings has been © 
«If the differential equation for the deflection curve of the beam on an elastic 
foundation is solved by the classical method, the number of simultaneous equa-— 

tions to be dealt with would be 4 n, where n is the number of discontinuities in 
loading. This, as the author himself has pointed out, is quite cumbersome, bel / 
is rather unfortunate that Penzienhas not mentioned in his introductionthe ex- 
istence and utility o of the operational methods in solving such discontinuity prob- — 
lems, With this method, the number of simultaneous equations to be dealt with Lo 
_ is only two, whatever the number of discontinuities, A paper onthe application — a 

of the operational method to simple beam 5 problems by William T, Thompson is r 
available.> Recently the writer® applied the operational technique to problems — 


of beams on elastic foundations in 1 which b beams with h arbitrary discontinuous 
_ By applying the Laplace transformation, the. ordinary differential ‘equation 
is reduced toa simpler form, usually a linear algebraic equation. This simple i. 
~ equation is then solved for a transition function from which the desired solu- 


“tee to the original differential equation is obtained by the use of inverse Lap- > 


lacetransforms, 


In addition the load functions that are used in the transformation method > 
provide, in the case of beams with arbitrary discontinuous loading, a very con- 
venient tool for arriving at the solution, in a continuous form, and in a single 
process, for moment, and reaction at any point ‘in the 


_ There is one more point of particular interest and that is va contrast to the 
4 n simultaneous eqeations, involved in the classical method, which stereo mgd 


_ tions, whatever the number of 


4 April, 1960, by Joseph Penzien, 
4 Research Fellow, Civ. and Hydr. Engrg. Sect., Indian Inst. of Scienc 
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Discussion By G Gordon K. Ray, F, ASCE, and William G. Westall, 9 


- advocating the use of steel reinforcement in ‘concrete airport pavements, ‘te 
we te 


: author draws upon three principal sources for his meee _ These are: 
Highway pavement performance. 


a * Corps of Engineers tests on reinforced concrete pavement. a 
7 
The purpose of this discussion is to ) objectively examine each of these, | 
keep the subject in proper perspective, it must be recognized that the use 
4 of reinforcement in highway pavements is but one of several methods being tried — 
by highway engineers it int their continuing s search for a + Gonign that will provide 


well stated by an ACI committee in its report on concrete slab dimensions:10 


9 a consideration of several factors, some eof which depend upon lon-— my 
gitudinal slab stresses, however produced, and others uponthe pavement = e 
. ee construction, maintenance, and past performance u under traffic, 
All of the factors involved must be so balanced against each other that 
the slab length finally chosen will result in optimum pavement perform- r 
ances, The economic phase of the problem is important and optimum 


performance must be attained at a cost which, although not necessarily 


The fact that several state highway departments build reinforced concrete 

in relatively long panels does not necessarily validate this practice under the 
criteria quoted previously. Using cement-treated subbases, California and 

several other states are building concrete pavements with relatively short joint 

_ spacing and without dowels for load transfer. The first pavement built on this 

_ design in California has shown excellent performance for 14 yr.*++ Neverthe- 

— less, the statement used previously also applies here: the successful use of 


this design does not validate it as the optimum, It must be assumed that each | 
of the many design methods now used ™ highway construction agencies was 


Mgr., Paving Bureau, Portland Cement Assoc., Chicago, Ill, 
9 Supervising Engr., , Airfield Program, Portland Cement Assoc, , Chicago, ni, 
“ 10 “Considerations in the Selection of Slab Dimensions,” ACI 53-24, ‘American Con- 
_ * il “Construction Practices on Cement-Treated Subgrades | for Concrete Pavements,” 
eem, presented at 39th Annual Meeting, Research Bd, 
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 Sanuary, 1961 
pare in an attempt to obtain 2 a balanced relationship between performance 
and economy. It appears likely that results from the AASHO Road Test near ~ 
Ottawa, Ill. will cause revision of some current highway design practices, + 
_ Following a comparison of joint spacing in reinforced and non- reinforced — 
pavements, the wine states in the last of his introduction: 


-@Both typ types s of design have been used and are now now being use used for air- 
‘port pavements, Reasons for using or not using reinforcement are the 
_ same as those for highways. Of course, they should be. The volume — 
_ traffic on airport pavements is much ~_— than on highways but the wheel 
loads are many ti times greater.” 


The in are subject to question, Because 


of a combination of factors, the reasons for using reinforcement in highway — 
pavements cannot be applied without exception to airport pavements, — of — 
these factors will be discussed because of their s special significance. Ly o ps 
Compton has pointed out that the traffic volume is much lesson airport — 
~ ments than on highways. ‘This difference in in frequency of loading will bear closer 


Airport traffic controllers have estimated 1 that, even under ideal conditions, - 


1, 000 to 1, 200 operations per day would — the ‘maximum rate of utilization at 


because it is based on the highest vaso rate possible, It does not take into 
account the reduced rate during weather minimums, nor allow for the slack 

hours of operation, In actual practice, the number of daily operations will be ‘ 

considerably below the theoretical maximum. Even if it is assumed that an 


airport pavement might be trafficked 1,000 times daily, this is less thanthe _ 


hourly rate on a busy urban. expressway or freeway. © _ Whereas it is true that 
only a relatively small percentage of thetotal number of vehicles on the high- 


is: also true of the aircraft thatuse airport pavements, Considering their much | 


way will equal or exceed the loading for which the pavement was designed, this — : 


_ greater traffic volume, highway pavements probably are subjected tothe — 
- mum allowable loads with greater frequency than airport pavements are sub- 
cow The manner in which loads are applied to pavements is a critical factor in 7 
the development of designcriteria, Westergaardand others have demonstrated | 
that loading at the free edge of aconcrete pavement creates stresses of greater — 
magnitude than interior loading where pavement continuity is provided by load 
transfer between joints. Airport pavements are seldom, if ever, subjected to 
Bonar along a freeor unsupported edge, On highways, however, wheels of the 
aviest vehicles travel close to the outer edge of the pavement, aos yea « 


‘ment edge, it must be concluded that highway pavements are subjected to more © 


severe conditions of loading than are airport pavements, The reasons for using 

The refers frequently (and logically) to Corps “ test data, 
: During the past twenty years the Corps has invested more than $10,000,000 in 
research and engineering investigations for the development of valid | design - 
‘procedures | | for airfield pavements, During the same period of time, as the 
construction agency for the Air Force, they have designed and constructed some 
$3, 000, 000, 000 worth of military airfield pavements, Probably the the — of 


he 
Considering their high traffic volume with load channelization at 
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‘See has 1 more experiencethan any other agency in the design, construc-_ 


tion, and performance evaluation of pavements built for modern aircraft. a Salas, 
is interesting, then, to look at their reasons for performing tests on 
-inforced concrete pavements, Ina paper . James P, Sale and — L. 


in “The increased thicknesses (up to inches) for the 100,0 
4 ound loading at the higher number of coverages being applied, 


"nesses s through steel or prestressing of the pavements. 
, Although it was realized that generally a reinforced rigid pavement — 
would be more costly than an equivalent strength non-reinforced rigid - 
Pav ‘ement, there was a definite need for such a design criteriato a 
some way in which certain pavement thicknesses could be reduced to meet 
surface grade requirements for different thicknesses in the various traf-_ 
and, when strengthening an existing pavement, it was impossible 
ws meet surface grade requirements without either removing and Teplac-— a>. 
ing areas of pavements or or over- -designing c certain areas. 


he 
Significantly, the ‘reinforced design cz criteria was developed for the purpose 
_ of meeting certain special conditions and no unsatisfactory performance of. 
plain concrete pavement was cited, On the contrary, this statement was in- _ 


cluded inthe same report: Ries 


quacy fromthe of construction feasibility, structural capacity, 
_ surface smoothness and durability for support of such propeller- -driven — 
aircraft as the B-29 and B-36 and more recently for the B-47 and -” 52 nae 


In view of their vast experience in all phases of engineering relating to air- 
_ field pavements, , it is of further interest to note that the Corps of Engineers E 
design manuall3 (with changes to March 1960) requires reinforced concrete — 
pavement only to meet certain conditions and not as general design criteria, 2 
ry Much of the author’s argument is predicted on assumptions regarding the be 
behavior of concrete pavements built without reinforcement or dowels. —iItis— 
implied that these pavements can be expected to develop uncontrolled cracking, 
- joint faulting, and general deterioration, This behavior pattern has not been _ 
confirmed by engineers who have made performance evaluation surveys of air- 
port pavements that areinservicen 
Plain concrete pavements that were adequately designed and built under good — 
construction procedures have not developed harmful cracks even under pro- 
Berns usage at capacity loading. This is well illustrated in the Corps of En- - 
t 14 on the enn, of test sections at Sharonville, Ohio, _ ; 


* Sale and Ronald L. Hutchinson, Journal of of the Air Transport Div Div. ASCE, 


13 “Engineering and Design, Rigid Airfield Pavements, ” EM 1110-45-303 


Engrs., , Dept, of the Army, February,1958, = 
an “Heavy Wheel Load Traffic on Concrete Airfield Pavements,” by Frank M. Mell- 7 
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These te test sections, which included various thicknesses of a and 

= _ plain concrete pavements, were subjected to 90,000 repetitions of a 100,000-Ib 
twin-wheel loading. The under-designed sections of each type of pavement 

= cracked and broke up under the load repetitions. The fact that the en 

* ‘pavement survived longer after the appearance of a crack seems" to have led on 

x to conclusions regarding the benefits of reinforcement. However, two of the —- 
non-reinforced sections were the only test pavements which did not crack under _ 
a number of load repetitions considered the equivalent of 10 yr of —o 

i we It has been demonstrated that dowels across contraction joints improve the 

of highway pavements that carry a large volume of heavy 


_ Regarding the assumption that dowels are required to prevent joint faulting 

in concrete airport pavements, reference i: is again made t to Corps of Engineers 

- design criteria which reflects the information collected in their long-term 
study of pavement performance. Th a The aaa design manuall3 gets forth these _ 


“Dowels w will be required across last dummy 
_ joints back from the ends of all runways, andsimilar dowel requirements — 
may be included for the ends of other large paved areas where local ex- 


perience indicates that they are needed.” 
It is true that the Federal Aviation Agency, requires dowels across con-_ 


traction joints in aprons, taxiways and the ends of runways “to provide an n- 


creased marginof safety with respect to load transfer in these critical areas.” a 
This appears to be a carry-over from their 1948criteria which was published 7 
at a time when the effect on pavement of multiple-wheel landing gears was not 
well understood. The cost of installing dowels, except at free pavement ends, - 
_ cannot be justified based on airport pavement performance. Within the past 
several years, the writers have made detailed surveys of millions of square — 
; yards of concrete airfield pavements designed and constructed by the Corps of 
Engineers, These pavements supported all types of military aircraft including» 
_ the heaviest jet bombers. Even on Air Force training | bases, heavy 
bomber traffic was at the saturation p point, there was no evidence of contraction 
_ joint faulting or other pavement defects that could be attributed tothe omission 4 
of dowels across contraction joints, Built without expansion joints except at 
abutting structures, airport. pavements a are restrained from movement to’ such 
; F an extent that aggregate interlock is retained at contraction joints. = = © 
Basedon his assumption that contraction joints inairport pavements 


dowels for proper the author ane the hypothesis that a 


of airport pavements built without. dowels at the, contracti ion joints, It was stated, 
also, that a saving in maintenance costs would result from reinforced pavement 
with fewer joints to maintain, +The fallacy of this theory has been conclusively | 
‘demonstrated by the performance of joints placed between long reinforced © 
panels. _ Excessive movement, ‘resulting in large joint openings, has made it 
to keept the joints sealed. TI The ‘infiltration of noncompressible ma- 


Aeronautics Administration, 
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terial into the wide joint openings has s caused spalling of concrete ; along the 

edges, Attempts to maintain the joints under these conditions have resulted in 

more maintenance effort than is required for | the normal joint spacing in non- 

. The plain concrete pavements built under Corps of Engineers criteria have 

5 _ provided the clean, smooth surfaces required for the operation of modern jet 
aircraft. Joint openings are minimized by the shorter panels and the infiltra- 7 

= tion of foreign materialis prevented, The ound of jet engine damage from > 

Spalled particles is thus eliminated, 


The great weights and peculiar g gear configurations of aircraft have 


caused the Corps of Engineers to become concerned over the possibility that 
_ the pavement thicknesses ~ required may eventually exceed the | limits of con- 
struction feasibility. This c concern need not extend to the pavements required 
for civil aircraft, ‘The current jet transports, oqnigges with dual-tandem landing ~ 


gears, require only about > 10 in, to 14 in, of plain concrete pavement even in— 
the critical areas of | loading su such as -taxiways ; and runway ends, In pavement | a 
strength criteria for the various classes of airports, the FAA has seenfit to 
provide for pavement strength considerably above the current aircraft require- 
ments, This provides a form of insurance against future increases in aircraft Sf 
a For example, at Dulles International Airport, the critical areas of 
concrete pavement are 15 in, in thickness, This provides sufficient pavement 
strength to support operations of a 500, 000-1b aircraft equipped with the dual- 
tandem type gear of the 300,000-1lb current jet transports. The possibility of 7 
these pavements becoming overloaded within the foreseeable future is extremely 
remote, The pavement thickness is well within the construction capabilities of 
_ This discussion is not intended in any way to disparage the value of steel 
reinforcement in concrete pavement where conditions justify its use, There 
is little evidence that the additional cost of reinforcement can be ye justified for : 
airport pavements, The use of long reinforced panels may be in some ways. 


objectionable from an operational standpoint. Properly designed plain concrete 
pavement has been proved teed adequate by many years of performance, 2 
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PROBLEMS OF STEEL MEMBER J 


0. ASPLUND, 15 ASCE. — _The Engesser- Vianello has proved 
quite useful in solving elastic buckling p problems, The work | needed depends o: on 
the number of iterations required and that is determined by the precisionof the | 
terminating estimate (Eq. 29(b) ) of the buckling load, This estimate is of gen-— 
eral importance and it may perhaps be ‘appropriate eo: discuss it herein, a 
a Generally speaking the buckling deflections n of an arch are soluble from a 


= 
g(n) + +Hh(n) =0. 
_ with boundary conditions; g and h are self- adjoint linear differential = 


ions with variable coefficients, For certain “ “eigen- -values” of the 


tal force H this boundary value problem | _ non- zero solutions . ie 


_ In the Engesser- Vianello method a roughly chosen approximation n stion (1) we 
-isfying the boundary conditions is substituted into the second term of Eq. 56 
"and an improved solution n= is by numerical or graphical inte- 


e(n(2)) = 
and the conditions. This new function n 2) | is asa ap- 
proximation and again substituted in the second term of Eq. 56. The solution © 
‘AL A. ated the convergence of such iterations, 
_ ing ¢ a: J. J. Koch18 proved the intuitive formula 


‘The | “eigenfunctions” ny are orthogonal and ca can be normatinet: a 


i nj) dx aes (60a) 

- 


June, 1960 by Georg Wastlund, 1 

15 prof,, Chalmers Inst. of Tech., Gothenberg 
F, Engesser, Ost. Ing. u. Arch, Ver, 1893, p. 506. 
pe A, Schwartz, Ges, Math, Abh, I, Berlin 1890, P. 241-265, 
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where = for j one 5 ij =j. All integrals here written are 
to be extended from x = 0 to x = 
Any roughly posed displacement nl) the boundary 
can now be expressed by a Fourier series 


Substituting 61 Land 57 57 into Eq. 58 


or if 7‘ (2) satisfies the contions 


Repeating this we 
= n 
and 
n-1 A(t) 


When n is large and M that is the case for fixed arches, 


> term of Eq. 65 needs to be considered and we find, proving Eq. 59, 
a Still, after a limited number nof iterations, the two. curves (n-1) and n(n) : 
are not perfectly affine, so the value (Eq. 59) of the eg safety factor will 
hich ent and 7\"/ are 
measured, It has therefore been suggested to use the ratio of the integrals of 


(n 


(n) =n (n- 
To avoid this, absolute va values of 7 can be used in Eq. 67. However, to _ cer-— 
tain that a correct limit is approached the writer suggests adopting, for ter- 


Inser erting the first ‘Eq. 65 we 


(1 ) 


4 


tions of with known exact and found that it is very easy to 
apply numerically and that it yields surprisingly accurate results, These ex- 
amples show, for instance, that two iterations using ‘Eq. 69 gives a more ac- 
curate estimate than three iterations using Eq. 


— 
— 
Gt ged s the best. The lowest buck- 
ee 0 and by Eq. 65 we obtaina a, 
— 
— 
-\2n-2. (1) Sn? im 
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‘= DYNAMIC RESPONSE OF ELASTU PLASTIC FRAMES@ 


AA. EREMIN, * M. ASCE.- Penzienstated, “Stiffness pro- 
_ vided by ‘secondary ‘components such as walls, partitions, stairwells, and so 
on, are primarily responsible for the apparent —— of many existing _—s 


_ tures to withstand strong motion earthquakes,” 
- It should be noted that when rigid walls or cutee are used as = 


members in a structure they should be designed to give a uniform variation of 


elasto-plastic rigidity in the structure, In Fig, 2 there has been shownthe dia- 


‘gram of elasto-plastic stress-strain variation, It would be interesting if rll 


would indicate the repeated stress-strain loop that resulted form a oa 
Generally it is considered that the strain in steel beyond the 


vy 


vibrating force, 
yield point does not recover in a short period of time, imi 


The graphsof dynamic stresses superimposed over the spectra of. stresses. 


for the El Centro earthquake in Fig. 5 to 12 isa biemacenagd contribution, It is 


~ Assoc. Bridge., Br — Dept., 
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119, In line 2 of Eq. change . Inline 3 of Eq.41 


121. line 2 of Eq. 47 change to In line 2 of Eq. 48 remove re the 
- first parentheses and insert a parentheses prior to the term 1/A. In line 5 of | : 
Eq. 48 insert a bracket just after the first minus sign and another bracket just 
ot to the first bracket that w er originally. In line 5 of Eq. 48 change 

the denominator nz E Ly ton E Iky . In line 6 of Eq. 48 insert a oe 1 


08 — the first plus sign and and ‘insert a bracket just prior to the the term ‘Sin ‘. 


1960 
—" 111. Delete the first paragraph o of the discussion by Mr. Newson and re-- 


_ Mr. Bonasia has presented avery interesting paper on a subject which 
is current and has especial significance for the younger engineer because 
it is fairly easy to understand; and, if called upon to do so, the younger 


engineer should have no dates solving the kind of — with which 


October, 1 


57. 7 should be written 


Eq. 35b should 


2k, Ke bs + (Kea Ate cd +4 + + Ket) 8 c+ 2 Ke f 


‘Eq. 35d should be Eq. a. 35¢; Eq. should | be Ea. 35d. 5d. 
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69-70. Eq. 36 ‘should be written = follows: 


| 
ab ved 


102. Change “react “reaction” to section” ion” 


p. 69. The proper Fig. 8 should be as follows: q 
— 
73, Eq. 43 should be writtenas | 
p. 88. Change Eqs. la, 1b and Ic to, = 
i dA=tax | (a) 
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-(HYS), 2479(SA3), '2480(STS), 2481(SA3), 2482(CO2), 2483(CO2), 2484(HY5), 2485(HYS), 2486(AT1)C. 2487 
1)¢, 2488(CO2)¢, 2489(HY5)©, 2490(SA3)°, 2491(STS)©, 2492(CP1), 2493(CO2),§ 
JUNE: 2494(IR2), 2495(IR2), 2496(ST6), 2497(EM3), 2498(EM3), 2499(EM3), 2500(EM3), 2501(SM3), 2502 
 (EM3), 2503(PO3), 2504(WW2), 2505(EM3), 2506(HY6), 2507(WW2), 2508(PO3), 2509(ST6), 2510(EM3),2511 

(EM3), 2512(ST6), 2513(HW2), 2514(HY6), 2515(PO3), 2516(EM3), 2517(WW2), 2518(WW2), 2519(EM3),2520 

(PO3), 2521(HY6), 2522(SM3), 2523(ST6), 2524(HY6), 2525(HY6), 2526(HY6), 2527(IR2), 2528(ST6), 2529 

JULY: 2541(ST7), 2542(ST7), 2543(SA4), 2544(ST7), 2545(ST7), 

2650(SU2), 2551(HY7), 2552(ST7), 2553(SU2), 2554(SA4), 2555(ST7), 2556(SA4), 2557(SA4), 2558(SA4), 
‘AUGUST: 2564(SM4), 2565(EM4), 2566(ST8), 2567(EM4), 2568(PO4), 2569(PO4), 2570(HY8), 2571(EM4), 
2573(EM4), 2574(SM4), 2575(EM4), 2576(EM4), 2577(HY8), 2578(EM4), 2579( PO4), 2580 
 (EM4), 2581(ST8), 2582(ST8), 2583(EM4)¢, 2584(PO4)¢, 2585(ST8)¢, 2586(sM4)¢, 2587(HY8)¢, 
SEPTEMBER: -2588(IR3) , 2589(IR3), 2590(WW3), 2591(IR3), 2592(HW3), 2593(IR3), 2594(IR3), 2595(IR3), 2596 
(HWS), 2597(WW3), 2598(IR3), 2599(WW3), 2600(WW3), 2601(WW3), 2602(WW3), 2603(WW3), 2604(HW3), 

2605(SA5), 2606(WW3), 2607(SA5), 200 KSA5)°, 2610(IR3), 2613(IR3)°, 

2615(EM5), 2016(EMs), 2617(sT10), 2618(SM5), 2620(EMS5), 2621(ST10), 

-:2622(EM5), 2623(SM5), 2624(EM5), 2625(SM5), 2626(SM5), 2627(EM5), 2628(EM5), 2629(ST10), 

9630(sT10), 2631( POS)¢, 2632(EM5)°, 2633(ST10), 2634(ST10), 2635(ST10)¢, 2636(SM5)¢. 
4 NOVEMBER: 2637(ST11), 2638(ST11), 2639(CO3), 2640(ST11), 2641(SA6), 2642(WW4), 2643(ST11), 2644(HY9), 

:2645(ST11), 2646(EY9), 2647(WW4), 2648(WW4), 2649(WW4), 2650(ST11), 2651(CO3), 2652(HY9), 2653(HY9), 


(EM6), 2676 (IR4), 2677 (HW4), 2678 (ST12), 2679(EM6), 2680(ST12), 2681 (SMB), 2682 (IR4), 
 2684(SMB), 2685(IR4), 2686(EM6), 2687(EM6), 2688(EM6), 2689(EM6), 2690(EMB), 2 2692 


2699(PP1), 2700(HY1), 2701(SA1), 2702(SU1), 270%(ST1), 
2707(HY1), 2708(HY1), 2709(POl1), 2710(HY1), 2711(HY1), 2712(ST1), ’2713(HY1), 2714(PO1), 2715. 
(ST 1), 2716(HY1), 2717(SA1), 2718(SA1), 2719(SU1)°, 2720(8A1)°, 2721(8T1), 2724 
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